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Abstract 

In a recent paper, we suggested that the density fluctuation spectra arising 
from power-law (or extended) inflation, which are tilted with respect to the 
Harrison-Zel'dovich spectrum, may provide an explanation for the excess large 
scale clustering seen in galaxy surveys such as the APM survey. In the light 
of the new results from COBE, we examine in detail here cold dark matter 
cosmogonies based on inflationary models predicting power-law spectra. Along 
with power-law and extended inflation, this class includes natural inflation. 
The latter is of interest because, unlike the first two, it produces a power-law 
spectrum without significant gravitational wave production. We examine a 
range of phenomena, including large angle microwave background fluctuations, 
clustering in the galaxy distribution, bulk peculiar velocity flows, the formation 
of high redshift quasars and the epoch of structure formation. 

Of the three models, only natural inflation seems capable of explaining the 
large scale clustering of optical galaxies. Such a model, though at best marginal 
even at present, has some advantages over standard CDM and on most grounds 
appears to perform at least as well. Power-law inflation's primary interest may 
ultimately only be in permitting a larger bias parameter than standard CDM; it 
appears unable to explain excess clustering. Most models of extended inflation 
are ruled out at a high confidence level. 
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1 Introduction 



Although recently it has been running into trouble from several sources, the standard 
Cold Dark Matter (CDM) model (White et al 1987; Frenk et al 1988; Efstathiou 
1990) has shown some considerable success in explaining the clustering of galaxies 
we see around us. This model is based upon the idea that there exists a primor- 
dial spectrum of density fluctuations which takes on a particular form known as the 
Harrison-Zcl'dovich or flat spectrum. Most frequently, this choice is motivated as 
being a prediction of the inflationary universe model. Nevertheless, it has long been 
known, at least within the inflation community, that the predicted spectrum is only 
approximately flat, and that within any given model there are readily calculable de- 
viations from scale invariance. An example is provided by the standard 'chaotic' 
inflation scenario (Linde 1983, 1987, 1990), wherein a logarithmic correction (see eg 
Salopek, Bond & Bardeen 1989; Mukhanov, Feldman & Brandenberger 1992; Schae- 
fer & Shafl 1992) adds some power at large scales, giving a density contrast at horizon 
crossing of 5-10% greater at lOOO/i"-*^ Mpc than at lQh~^ Mpc. Hitherto, when our 
knowledge of the primeval spectrum was restricted to studies of local clustering up 
to scales of perhaps 100/i~^ Mpc, these corrections where rightly regarded as insignif- 
icant. The recent results (Smoot et al 1992) from the Cosmic Background Explorer 
satellite (COBE) have changed this picture, by providing for the flrst time good es- 
timates of the amplitude of the spectrum on scales of upwards of 1000/i~^ Mpc. The 
time has come for these correction terms to be taken seriously, and possibly the real 
test of inflation lies not in how flat the primeval spectrum is, but in the type and size 
of deviations from it. 

A certain class of inflationary models predict an even more dramatic change 
from the Harrison-Zel'dovich case, wherein the primeval spectrum generated exhibits 
power-law deviations from flatness. There are several such inflationary scenarios. The 
simplest is power-law inflation (Abbott & Wise 1984b; Lucchin & Matarrese 1985), 
which can be realised via a scalar field evolving in a potential of exponential form 
and which leads to an expansion of the universe which is not the conventional nearly 
exponential growth, but rather has the scale factor growing as a rapid power law in 
time. A second option is provided by the so-called extended inflation scenarios (La 
& Steinhardt 1989; Kolb, Salopek & Turner 1990; Kolb 1991), wherein modifications 
to Einsteinian gravity allow infiation to proceed via a first-order phase transition. 
Yet another option is provided by 'natural inflation' (Freese, Frieman & Olinto 1990; 
Adams et al 1992), a speciflc case of the more general situation of a scalar held 
evolving near the top of an inverted harmonic oscillator potential which also gives a 
power-law spectrum. 

While power-law primeval spectra are often written down, the consequences of 
such a spectrum in the complete CDM cosmogony have not been fully investigated. 
The first detailed discussion of such models was made by Vittorio, Matarrese and 
Lucchin (1988), where they also allowed the density parameter Q, to be less than one. 
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Subsequent investigations by these authors and collaborators (Tormen, Lucchin & 
Matarrese 1992; Tormen et al 1992) have concentrated primarily on analysis of bulk 
velocity flows. Interesting constraints in the n-ft plane were also obtained by Suto 
and Fujita (1990) and by Suto, Gouda and Sugiyama (1990) based on predictions for 
the cosmic Mach number combined with the then current microwave anisotropy lim- 
its. A pre-COBE analysis much less specific than that given here has been given by 
Kashlinsky (1992). In a recent paper with Will Sutherland (Liddle, Lyth & Suther- 
land 1992, henceforth LLS), we examined the consequences of such spectra for galaxy 
clustering. The motivation lay in recent observations such as the APM survey (Mad- 
dox et al 1990, 1991), which have shown that on large scales there appears to be 
stronger clustering than the standard CDM model can provide. These power-law 
spectra do indeed possess extra power on large scales. We found that for a certain 
range of values of the slope, one could reproduce the clustering seen in the APM 
survey, while withstanding the then current microwave background bounds. 

Given this success, the emergence of the COBE data makes this an excellent time 
to examine this scenario in more detail, especially as in a certain sense the COBE 
data too show a rather larger amplitude of fluctuations than one would expect in the 
standard biassed CDM model. We have also taken this opportunity to analyse several 
other aspects of the scenario, though we have not gone so far as to carry out iV-body 
simulations. Instead, we have restricted ourselves to calculations which can be done 
analytically or semi-analytically. 

Recently there have been several papers analysing various aspects of power-law 
spectra (Salopek 1992; Davis et al 1992; Liddle & Lyth 1992; Lucchin, Matarrese 
& Mollerach 1992; Adams et al 1992; Lidscy & Coles 1992), which have however 
tended to deal with specific issues relating to particular models. With the exception 
of Adams et al, these papers have also failed to come up with detailed constraints on 
the scenario. Our aim here is to provide a more unified treatment, covering all the 
infiationary possibilities for generating power-law spectra, with the ultimate aim of 
constraining the slope of the spectrum and the bias parameter for the different mod- 
els. One has already seen limits on the slope quoted by COBE (Smoot et al 1992) of 
n = 1.1 ± 0.5, these numbers being simply derived from the variation of microwave 
anisotropics with scale between around 10° and the quadrupole, and thus being inde- 
pendent of the specific choice of dark matter. [It has been suggested (D. T. Wilkinson, 
conference talk, Oxford, June 1992) that a more detailed error analysis (recognising 
that the error bars on the data are not independent) will increase the size of the 
errors on n.] It is natural for one to expect that if one makes the additional require- 
ment, within a specific cosmogony, that the appropriate galaxy clustering is obtained 
at scales a factor of 100 smaller than those COBE samples, then the hmits on the 
spectral slope become yet tighter. On the other hand, one should bear in mind that 
the COBE error bars on n arc presumably, like the rest of the data, 1-sigma error 
bars and so a stronger level of exclusion would allow further deviations from scale 
invariance. 
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Our general aim is to constrain the spectral index and bias parameter in the 
cold dark matter cosmogony by finding limits on the spectral amplitude at difi'erent 
length scales. Of most interest for this purpose will be the amplitude on large scales 
(> 1000/i~^ Mpc) from the COBE experiment, and on intermediate scales using the 
velocity data (which see the mass spectrum itself rather than the galaxy spectrum). In 
particular we shall use results from the QDOT survey (Saunders et al 1991; Saunders, 
Rowan-Robinson & Lawrence 1992) which fix the amplitude at around 20/i~^ Mpc. 
It is important to note that both these scales are expected to be well into the linear 
regime and thus linear calculational techniques can be used. To aid comparison with 
other work, we shall normally give the inferred normalisation in terms of the optical 
galaxy bias parameter h defined later, which is assumed to be the reciprocal of the 
mass dispersion in 8/i~^ Mpc spheres using the CDM spectrum. Note that this use 
should not be taken to imply that the results have anything to do with optical galaxies 
or scales which are close to entering the nonlinear regime. 

Having obtained limits in this way on the mass spectrum, reference to optical 
galaxies provides limits on the bias parameter, enabling us to present our constraints 
in the form of allowed regions in the n-h plane. In order to clarify why particular 
choices of n appear in these sections, let us presage here the ultimate result that n 
is constrained in any infiation model to be no less than 0.7, at very high confidence. 
With the allowed values in mind, we then go on to investigate whether any further 
constraints can be imposed by an examination of aspects of the formation of non- 
linear structure, wherein however the calculations are technically more tricky and 
comparison with observation much more ambiguous. 

The layout of this paper is thus as follows, and generally works from large to small 
scales. Firstly, section |^ discusses the infiationary models giving rise to power-law 
spectra. In section H, we discuss our processed CDM spectrum and its normalisation. 
Our aim here is to be as explicit as possible in our definitions and in where the 
uncertainties lie. Section ^ discusses large angle microwave background anisotropies, 
a discussion dominated of course by the observations of the COBE satellite. These 
calculations refer to the largest scales inside our horizon. In section ^, we review the 
results of LLS concerning the application of power-law spectra to the excess large- 
scale galaxy clustering. In this section we also discuss the bulk fiows, which serve as 
another indicator of the power spectrum amplitude on these scales. Section |^ deals 
with the issue of shorter scales, corresponding to typical galaxy and cluster masses. 
There we investigate a range of topics concerning the entry into the nonlinear regime 
In this region of the spectrum, our power-law spectra offer less power than standard 
CDM, and hence one expects the structure formation constraint to be more stringent 
here. Finally, we conclude and discuss our results. 

The universe is taken to have critical density, with present Hubble parameter 
Hq = lOO/ikms"^ Mpc where h = 0.5. For ease of comparison with other work we 
often take Mpc for the distance unit. The units are taken to be such that c = 1 
and (in section ||) = 1, and the Planck scale is defined by mp; = 
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2 Inflationary Models giving Power-Law Spectra 




As discussed in the introduction, there are at least three different inflationary models 
giving rise to power-law spectra. While there are certainly relations between these as 
we shall see, each has its own features and we shall discuss each in turn. 

Let us first briefly sketch the derivation of the spectra, following Liddle and Lyth 
(1992). A scalar field (f) moving in a potential V is usually taken to satisfy slow- roll 
conditions (' = d/d0) 

(2.1) 

^ (2.2) 

H = a/ a is the Hubble parameter during inflation, where a is the scale factor nor- 
malised to a = 1 at present. Denoting the epoch at the end of inflation by a subscript 
2, and the epoch k = aH when a scale k leaves the horizon with a subscript k, the 
number of Hubble times N{tk) between a scale k leaving the horizon and the end of 
inflation is 

In the slow-roll approximation, the spectrum of density fluctuations is given by 



PI 



where 5h is roughly the density contrast at horizon entry and is defined precisely in 
the following section. The efTective spectral index is given by 

1 - n ^ = (66, - 2r^,) (2.6) 

We see immediately that in order to have significant deviations from the fiat n = 1 
spectrum one must come close to violating at least one of the slow-roll conditions. 
Corrections to these slow-roll results have been calculated exactly for power-law in- 
flation (Lyth & Stewart 1992a), and to flrst order in the departure from slow- roll for 
the general case (Lyth & Stewart 1992b). As discussed in our earlier paper (Liddle 
& Lyth 1992) they are not important in practice. 

The models giving rise to power-law spectra are as follows. 
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1. Power-Law Inflation. 

Power-law inflation a (x arises (Abbott & Wise 1984b; Lucchin & Matar- 
rese 1985; Barrow 1987; Liddle 1989) when a scalar field (f) evolves down an 
exponential potential of the form 



It is of particular interest to researchers in inflation because although the scalar 
field comes close to violating the slow-roll conditions, exact analytic solutions 
exist both for the dynamics of infiation and for the density perturbations gen- 
erated. The exact formula for the spectral index is (Lyth & Stewart 1992a) 

n=l-2/(p-l) (2.8) 

In the limit of large p the spectrum tends rapidly towards the fiat n — 1 spec- 
trum. 

2. Extended Inflation. 

Extended infiation (La & Steinhardt 1989; Kolb 1991) looks very different from 
power-law infiation, but is in fact related during infiation by a conformal trans- 
formation which allows the usual power-law infiation machinery to be utilised in 
its study. It is based on modifications to the gravitational sector of the theory, 
which allow a first-order infiationary phase transition to complete satisfactorily. 
The original model (La & Steinhardt 1989) was based on a Brans-Dicke theory 
with parameter cu, and although this proved insufficient to allow present day 
tests of general relativity to be satisfied, it has remained the paradigm around 
which more complicated working models are based. The spectral slope in such 
a model is given by (Kolb, Salopek & Turner 1990; Guth & Jain 1992; Lyth & 
Stewart 1992a) 

2oo-9 , , 

n = 2.9 

2uj-l ^ ^ 

with extended infiation being conformally equivalent to power-law infiation with 
2p = a; + 3/2. 

A crucial difference between power-law and extended inflation is that extended 
inflation suffers an additional constraint, as one must avoid the large bubbles 
generated as inflation ends from being so profuse as to unacceptably distort the 
microwave background. This constrains a; as a function of the infiaton energy 
scale M as (Liddle & Wands 1991; Liddle & Lyth 1992) 

a; < 20 + 0.71ogio(M/mpi) (2.10) 
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Bounding M using the microwave background limits on the fluctuation ampli- 
tude gives a; ^ 17, corresponding to n ^0.75. We have discussed the signifi- 
cance of this bound in an earlier paper (Liddle & Lyth 1992); the prognosis for 
the extended inflation model is not good, as we shall recap in this paper. 

3. Natural Inflation. 

The natural inflation model (Preese, Prieman & Ohnto 1990; Adams et al 1992) 
is based on a pseudo-Nambu-Goldstone boson evolving in a potential 

V(0) = A^(l±cos(</.//)) (2.11) 

where A and / are mass scales. Evolution near the top of this potential (as is 
required for sufficient inflation with parameter choices significantly tilting the 
spectrum), gives rise to a power- law spectrum. In fact, this is a realisation 
of the more general case of a scalar field evolving near the top of an inverted 
harmonic oscillator potential 

V{(j)) = Vo- imV^ (2.12) 



which gives a spectrum of slope 



2 2 

m m PI 



n = 1 - — -/^ (2.13) 
For natural inflation one has n = 1 — m^^JSTip. 

That natural inflation gives the same spectrum as power-law inflation may seem 
surprising in the light of claims that one can reconstruct inflationary potentials 
from a given spectrum (Hodges & Blumenthal 1990). In fact, these two models 
can be regarded as different regimes of an all-encompassing potential (which 
one can calculate in the manner of Hodges and Blumenthal) which is essentially 
1/ cosh^ with various factors thrown in. In the ~ region we have the 
inverted harmonic oscillator, while at large (j) we have the exponential region. 
In a sense, this is the unique potential from which all inflationary models giving 
power-law spectra arise. With this potential, the scalar field can roll from 
the top down to the exponential region, while in a slow-roll approximation 
generating an exact power-law spectrum. However, one must also note that as 
discussed above, the fact that the spectrum has been tilted implies that the slow- 
roll approximations are at best only just satisfied, so there will be corrections to 
the slow-roll spectrum in all regions of this potential. In the exponential region 
these corrections are known to affect only the amplitude and not the slope; one 
hopes for the same near the top of the potential. 

Each of the above models can produce a power-law spectrum of the appropriate 
amplitude. However, as discussed in many recent papers (Krauss and White 1992; 
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Salopek 1992; Davis et al 1992; Liddle & Lyth 1992; Lucchin, Matarrese & MoUerach 
1992; Souradeep & Sahni 1992), one must pay attention to the predicted amplitude 
of gravitational waves (this was noted for power-law inflation long ago by Fabbri, 
Lucchin and Matarrese (1986)). In regimes where the slow-roll approximation is in- 
valid or close to violation, one expects that these can become large, and the above 
work has confirmed that for power-law inflation (and by inference for extended infla- 
tion), the microwave contributions from gravitational waves may not only be large, 
but can dominate those from the scalar density fluctuations. To a good enough ap- 
proximation for our purposes (see Lucchin, Matarrese & Mollerach (1992) for more 
precise results), the ratio R of the tensor to the scalar contribution to the squared 
microwave multipoles is independent of the multipole and given by (Liddle & Lyth 



where the potential and its derivative are to be evaluated at the point corresponding 
to a given scale passing out of the horizon. In general R is scale-dependent, but for 
power-law and extended inflation one has the scale-independent ratio 



and thus for values of p less than around 12 (n < 0.82) the tensor modes dominate 
the observed microwave background anisotropies. 

On the other hand, for natural inflation the ratio R is extremely small as the rel- 
evant scales leave the horizon, and hence there are no gravitational wave corrections 
of note in that case. Thus natural inflation (and similar models) provide an inflation- 
ary mechanism for generating power-law spectra without the large gravitational wave 
accompaniment. We shall return to the role of gravitational waves in section | for 
their microwave background implications, but for now we shall turn to the spectrum 
of density fluctuations. 

3 The Spectrum 
3.1 Definitions 

The matter density contrast is 5(x, t) = (p(x) — p)/p, where p(x) is the matter 
density, p is its spatial average, t is time and x = (x^,x^,x^) are comoving Cartesian 
coordinates. There is some freedom in the definition of t and x, but this 'gauge 
freedom' is not significant during the matter dominated era which concerns us (Lyth 
& Stewart 1990). As long as they are small, 6 and related quantities evolve according 
to the linear equations of cosmological perturbation theory, with each Fourier mode 
evolving independently. At each comoving point during the matter dominated era, 6 
is proportional to the scale factor a, or equivalently to {aH)~'^. 



1992) 




(2.14) 



R ~ 12.4/p 



(2.15) 
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According to the inflationary models described in the last section the primeval 
density contrast is a Gaussian random field at very early times, and therefore as long 
as linear evolution holds. As such its stochastic properties are completely determined 

by its spectrum. Different definitions of the spectrum exist in the literature. Letting 
5k denote the Fourier coefficient of the density contrast in a box of coordinate volume 
V, with the normalisation 5k = J e^''^'^6{x)(fx, we work with the spectrum V 
defined by 

V^V{ky2n'){\S^\') (3.1) 

where the bracket denotes the average over a small region of /c-space. With this 
definition the dispersion a of the density contrast is given by 

^' = IX (3.2) 

In the literature, our V is denoted variously by Vp (Salopek, Bond and Bardeen 
1989), Pp (Lyth & Stewart 1990), (Kolb & Turner 1990), 5pk/p (Linde 1990) 
and dap/din k (Bond Sz Efstathiou 1991). The quantity usually denoted by P in the 
literature is equal to Nk~^V, where N is an author-dependent normalisation factor 
which is often left undefined. 

The standard assumption for the /c-dependence of the primordial spectrum is 
P (X k, corresponding toV oc k'^. This is variously referred to as a Harrison-Zel'dovich 
spectrum, a flat spectrum, or a scale- free spectrum, and it is generated by the vacuum 
fluctuation of the inflaton field if the slow-roll conditions are well satisfied. Here we 
are exploring a power- law primordial spectrum P oc /c", corresponding to V (x /c^"*"", 
with n < 1. 

During the matter dominated era, the spectrum may be written 

nk) = (^-^^' nk)5Uk) (3.3) 

Up to a gauge dependent numerical factor (Lyth & Stewart 1990), 5jf is the primeval 
spectrum at the epoch k = aH when the scale k leaves the horizon. (In the literature 
Snik) is often written as Sp/p\hor ) Its scale dependence is 6jj oc A;"^"*^, which is why 
the standard choice n = 1 is referred to as a flat or scale-invariant spectrum. 

The quantity T(k) is known as the transfer function. It measures distortions 
to the primeval spectrum generated once a given scale comes within the horizon. 
On those very large scales which enter the horizon long after matter domination 
{k~^ ^ 100 Mpc) it is equal to 1. To calculate it on smaller scales one has to know 
the matter content of the universe (see Efstathiou (1990) for a discussion covering 
various physical situations). For hot dark matter, free streaming erases perturbations 
on scales up to around the horizon scale at matter-radiation equality, but for cold 
dark matter the predominant effect is simply that the growth of modes which enter 
the horizon during radiation domination is suppressed. 
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Taking the dark matter to be cold, the standard picture of the early universe 
determines a practically unique transfer function once the baryon density Jl^ is fixed. 
There are however considerable discrepancies between the parametrisations which 

appear in the literature. For Qb = ^ there are two widely quoted parametrisations, 
one due to Bond and Efstathiou (1984) and the other to Bardeen et al (1986). We 
use the former, as given by Efstathiou (1990) 

T{k) = [l + (ak + {hkf + {ckf) (3.4) 

where 

a = 6.4(n/i2)"^Mpc ; 6 = 3 (Q/i^) Mpc ; c = 1.7 (fi/i^) Mpc ; i/ = 1.13 

(3.5) 

In figure 1 this parametrisation is compared with that of Bardeen et al (1986), and 
with two others. At scales around k = 27i/10h^^ Mpc^^, the discrepancies can be of 
order 10%. In particular, the parametrisation of Davis et al (1985), which has the 
fewest fitting parameters, appears rather less accurate than the others. 

Adams et al (1992) claim that Bardeen et aFs parametrisation is extremely ac- 
curate for Qb = 0, and that increasing to the value Qb — -06 favoured by 
nucleosynthesis decreases it by about 15% at A; ~ 27r/10/i^^ Mpc. This would bring 
it 5% below the parametrisation that we have adopted, forcing the normalisation up 
by roughly that amount. 



3.2 The filtered density contrast 

The density contrast 5(x) will evolve linearly as long as its dispersion satisfies cr <C 1, 
except in those rare regions where it becomes ^ 1 and gravitational collapse takes 
place. Following the usual practice, we assume when necessary that the linear evo- 
lution is at least roughly valid right up to the epoch a ~ 1. Soon after that epoch, 
a large fraction of the matter collapses into gravitationally bound objects, and lin- 
ear evolution becomes completely invalid. The epoch of non-linearity, defined as the 
epoch when a — 1, corresponds to a redshift z^i given by (1 -|- ^ni) = l/^O; where (Tq 
is the linearly evolved quantity evaluated at the present time. 

A 'filtered' density contrast 5(i?/,x), which has a smaller dispersion a{Rf), can 
be constructed by cutting ofi^ the Fourier expansion of 5(x) above some minimum 
wavenumber ~ ^/Rf, or equivalently by smearing it over a region with size — Rf- 
The filtered quantity will evolve hnearly until the later epoch o-{Rf) — 1. 

A precise definition of the filtered density contrast is made by means of a 'window 
function' W{r), which is equal to 1 at r = and which falls off rapidly beyond some 
radius Rf (Peebles 1980, Kolb & Turner 1990). The filtered density contrast is 

S{Rf, x) = y W{Rf, |x' - x|)(5(x')(iV (3.6) 



9 



and its spectrum is 

where 

and 



V{Rf,k)= W{Rf,k)/Vf V{k) 



W{Rf,k) 



Vf 



The filtered dispersion is 

a\Rf) 



e~'^-''W{Rf,r)£x 
W{Rf,r)d^x 



W{Rf,k)/Vf] V{k)- 



k 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



The quantity Vf is the volume 'enclosed' by the filter. It is convenient to define 
the associated mass M = poVf, where po is the mean comoving density which in an 
h = 0.5 universe is given by po = SH^/SnG = 6.94 x W^Mq Mpc"^ (M© being the 
solar mass). One normally uses M instead of Rf to specify the scale, writing 6{M, x) 
and cr(M). 

The two popular choices are the Gaussian filter 



W{Rf,'. 



exp(-rV2i?J) 



Vf = {2T:f'^R) 

W{Rf,k)/Vf = exp{-kRf) 

M = 4.36 X lO^'^h^iRf/lMpcfMQ 

and the top hat filter which smears uniformly over a sphere of radius Rf 

W{Rf,r) = 9{r-Rf) 



Vf = A-rcRyS 



W{Rf,k)/Vf = 3 



'sm{kRf) 



cos 



jkRf) ] 



{kRff J 



M = 1.16 X 10'^/i^(i?//lMpc)^Me 







(3.11) 
(3.12) 

(3.13) 
(3.14) 



(3.15) 
(3.16) 

(3.17) 

(3.18) 



The Gaussian filter is the most convenient for the theoretical calculations described 
in Section but the top hat filter is widely used to as a means of presenting data as 
in Section |. 



3.3 The conventional normalisation 

The final ingredient required to fix the CDM spectrum is the overall normalisation. 
In the past this has usually been done by comparing the linearly evolved theory with 
observation, for the galaxy correlation function on a scale of order 10h~^ Mpc. Two 
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standard prescriptions exist (Efstathiou 1990). The first requires that the dispersion 
(7p(r) of galaxy counts in spheres of radius r has the observed value unity at r = 
Sh~^ Mpc (Davis & Peebles 1983). The second requires that J3(r), the integral of the 
second moment of the galaxy correlation function up to distance r, has the observed 
value (Davis & Peebles 1983) 270 h"-^ Mpc'^ at r = lOh-^ Mpc. 

These normalisation schemes both refer to the statistics of the clustering of lumi- 
nous galaxies, not to the clustering of mass as described by the spectrum V. It is 
an important ingredient of the CDM cosmogony that these are not necessarily taken 
to be the same — so-called biassed galaxy formation (Bardeen et al 1986; Efstathiou 
1990). To be specific the mass density correlation function ^, related to the spectrum 
by 

X /""^ _ ,, , sin A;r d/c 
e(r) = /_ V{k)-—- (3.19) 



/o kr k 

is assumed to be related to the correlation function of luminous galaxies by 



(3.20) 



with a bias parameter b which is roughly independent of the scale r. Including the 
bias parameter, the normalisations implied by the two schemes are 



{r)/b' 



POO 

= 9 / Voik) 
Jo 



sin kr cos kr 



1 2 



(kr)^ (kry 
at r = 8h~^ Mpc 



d^ 
T 



(3.21) 



and 



{.73/bf 



sin kr 



— cos kr 



at r 



kr 

IQh-^ Mpc 



d^ 
T 



(3.22) 



In each case, Vo is the linearly evolved spectrum, evaluated at the present epoch. 

One sees that (7g{r)/b is the dispersion of the linearly evolved density contrast 
with top hat filtering on the scale r, whereas {2r^'^ J'^ir))^^"^ /b is its dispersion on the 
same scale with a different choice of filter. As noted earlier, the dispersion of a filtered 
quantity evolves linearly only as long as it is ^ 1, and this condition is only marginally 
satisfied by the quantities encountered in the two conventional normalisation schemes. 
The non-linear correction will be worse if 6 ~ 1 than if 6 ~ 2, and CoTichman 
and Carlberg (1992) have claimed on the basis of numerical simulations that it is 
indeed significant in the former case. Evaluating the integrals, one finds that the 
J3 normalisation is typically 10% lower than the normalisation using cTg, across the 
range of n in which we shall be interested. These are shown in figure 2. As remarked 
above, different choices of transfer function parametrisation can give a correction of 
around 10% as well. 
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In what follows we always adopt the ag normalisation scheme, in the sense that 
normalisation of the spectrum is specified by quoting the quantity h defined by 
Eq. ( p. 21 ). In words, is defined as the present value of the dispersion of the linearly 
evolved matter density contrast, with top hat filtering on the scale Rf = 8h~^ Mpc. 
Only occasionally do we invoke the much stronger assumptions, that the galaxy cor- 
relation function is given by C,gg{r) ~ }P'^{r) and that ^(r) is equal to the linearly 
evolved quantity on scales r ~ lO/i^^ Mpc. Our most significant results, summarised 
in figure 11 later in the paper, do not involve these assumptions. 

Figure 3 shows hV^ {k) for a selection of n values, where as always V is the 
linearly evolved quantity and the subscript zero indicates evaluation at the present 
epoch. The spectra for different n typically cross at a scale of around = 2n/32h~^ 
Mpc. Figure 4 shows the corresponding dispersion ba{M) for n = 1 and n = 0.7, for 
both Gaussian and top hat filtering. The mass scale M runs from M = lO^Mo, the 
Jeans mass at decoupling (Peebles 1980, Kolb & Turner 1990) to M = lO^^M©, the 
mass of the universe on the scale of the CfA survey. The comoving wavenumber 
runs over the corresponding range of filtering radius Rf. 

For n = 1 the spectrum V{k) is practically flat on small scales. For n < 1 the 
power on small scales is reduced, though even for n = not by as much as in the 
standard hot dark matter model (White, Frenk & Davis 1983) or with string-seeded 
hot dark matter (Albrecht & Stebbins 1992)). As a result the spectrum starts to fall 
on scales k~^ < .IMpc. The implication of this change for small scale structure is 
discussed in Section 0. 



4 Large Scales: Large Angle Microwave Background 
Anisotropies 

Measurements of the cosmic microwave background anisotropy at large angular scales 
provide the only evidence we possess as to the amplitude of the spectrum on scales 
> 100h~^ Mpc. Ignoring the possibility of reionisation, the surface of last scattering 
for the cosmic microwave background lies some 190/i~^ Mpc inside the particle horizon 
(Hogan, Kaiser & Rees 1982), whose present distance is 2H^^ = 6000/i~^ Mpc. This 
means that an angular scale 6q corresponds to a linear scale r = 1016h~^ Mpc. 

For angular scales ^ 1°, the linear scale is much bigger than the horizon size 
at last scattering. This means that the microwave anisotropy does not involve such 
awkward physical processes as Thomson scattering from moving electrons and the like, 
which vastly complicate matters when one reaches scales of a degree or less (Bond 
& Efstathiou 1987). It is in fact given by the well-known and understood Sachs- 
Wolfe effect (Sachs & Wolfe 1967), which corresponds simply to fluctuations in the 
gravitational potential across the microwave sky caused by the density fluctuations. 

If re-ionisation occurs the distance of the surface of last scattering will be multi- 
plied by some factor e < 1, and the Sachs- Wolfe formula will be valid only on scales 
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3> degrees. Re-ionisation is not thought to be hkely in the standard CDM model, 
due to the relatively late onset of nonlinear behaviour which is further exacerbated 
with power-law spectra. Prom now on the possibility of re-ionisation is ignored. 

The only data we consider here are those from the COBE experiment (Smoot et 
al 1992), which to date provides the only positive detection of anisotropics. It corre- 
sponds to an angular resolution of order 10°, where one may safely ignore the pos- 
sibility of re-ionisation. The corresponding hnear resolution is of order 10^h~^ Mpc. 
Experiments with smaller angular resolution have so far reported only upper limits, 
which are compatible with the COBE detection for the CDM model with n <1. 

Let us first review the pertinent results from COBE (Smoot et al 1992). The most 
useful and significant measurement made is that of the variance of the temperature 
fluctuations smoothed on a scale of 10°. The value quoted is (1.1±0.2) x 10~^. Like all 
the COBE results, the error limits quoted are 1-sigma. Of less theoretical importance 
is the measurement of the quadrupole anisotropy Q2, which is also susceptible to 
experimental contamination from the galaxy in particular. This is quoted as Q2 — 



COBE also have results on all intermediate scales, readily obtained by the appro- 
priate smoothing of their basic data set. These are the data used to calculate the 
1-sigma limits on the best fit primeval spectral slope n = 1.1 ± 0.5. Unfortunately, 
this intermediate data has not been given in the paper, and so consequently we must 
work solely with the information given above and the slope limits. 

Now we proceed to the formahsm, which follows closely that of Scaramella and 
Vittorio (1990). [See also Abbott & Wise 1984a; Bond & Efstathiou 1987; ScarameUa 
&: Vittorio 1988; Efstathiou 1990, 1991.] For studies at large angles, the most con- 
venient approach is the expansion of the temperature fiuctuation field in spherical 
harmonics. Dropping the dipole term which is dominated by our peculiar velocity 
relative to the comoving rest frame, one writes 



where 6 and are angles on the sky and x is the observer position. With gaussian 
statistics for the original density field, the coefficients af'(x) are gaussian distributed 
random variables of position, with zero mean and with a rotationally invariant vari- 
ance which depends only on /, 



where the angle brackets are averages over all observer positions. It is important 
to appreciate however that what one can actually measure is only the temperature 
anisotropics at a single observer point, our own. Thus in practice one only gets a 
single realisation of the probability distribution functions given above. This leads to 
statistical uncertainties, as an individual realisation may not reflect the properties of 



(5 ±1.5) X 10"^ 




1=2 m=-l 



(4.1) 




(4.2) 
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the ensemble averaged system. This effect, normally called the cosmic variance, is 
significant for observables which depend only on a limited number of the a^, and in 
particular for the quadrupole which only depends on the five (two of which are re- 
dundant rotational information). This is the fundamental reason why the quadrupole 
measurement is not particularly useful for constraining theories. On the other hand, 
the variance at 10° depends on a significant number of the a™ and is far less suscep- 
tible to statistical vagaries. 

The multipoles Qi as seen from a particular observer point are defined from the 
variances of the individual of* as 

Qi-ii: (4-3) 

m=—l 

In fact, the temperature autocorrelation function measured by a single observer 

where the angle brackets indicate an average over all directions separated by an angle 
a, is given exactly in terms of the Qf via 

CO 

C{a) = Y^Q^Pi{cosa) (4.5) 

1=2 

Note that our Qf have a different normahsation to that of Scaramella and Vittorio 
(1990). It is defined here so that \JqI is exactly the quantity COBE measures. 

Since the are the sums of squared gaussian random variables, they are described 
by chi-squared distributions with the appropriate number of degrees of freedom; for 
the l-th multipole there are 21 + 1 degrees of freedom, and in particular the quadrupole 
is described by a x| distribution. 

The expectation of Qf, again averaged over observer positions, is directly related 
to the variance of the aj", via 

(Qf) = ^ (2/ + (4.6) 

The values of the Hf are readily calculated. On these large angular scales, the temper- 
ature anisotropics from the density perturbations are dominated by the Sachs- Wolfe 
effect (Sachs & Wolfe 1967). The appropriate formula for this curvature contribution 
is (Peebles 1982) 

foo A h 



E^(scalar) = tt — jf {2k/ aH) 5jj (4.7) 

where we write 'scalar' to distinguish from other contributions to the variances. Here 
ji is the spherical Bessel function. With 5'^{k) cx A;"~^ this becomes 



(scalar) = - 
8 



^/(/ + l)r((3-n)/2)r(/ + (n-l)/2) 



2 ' 'r((4-n)/2)r(i + (5-n)/2) 



2^ + ^-5|,(i7o/2) (4.^ 
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As discussed in section there is in power-law and extended inflation models an 
additional contribution to the multipoles from gravitational waves. The contribution 
to the squared multipoles from gravitational waves is adequately given (for the range 
of / values we will use) from the scalar contribution above just by 

S?(grav) = — S? (scalar) (4.9) 
P 

so that 

15 — 13ri 

E;^ (total) = — (scalar) (4.10) 

In contrast, natural inflation has a negligible gravitational wave contribution, even 
though it can give a tilted spectrum just like the power-law case. 

4.1 The quadrupole anisotropy 

As we shall see below, the quadrupole is not particularly useful for constraining the 
models, and we shall ultimately drop it in favour of the 10° result. Nevertheless, it is 
worth examining the predictions of the model to see why we come to this conclusion. 

The above equations readily allow us to calculate the mean of the predicted distri- 
bution function for Q\. It is given in figure 5, again as a function of spectral slope, for 
both types of inflationary model. Because the spectral normalisation is uncertain by 
the bias factor, we have plotted the quantity b'^{Ql). The COBE result of 2.5 x 10~^^ 
corresponds almost exactly to the mean prediction for the flat spectrum when the 
bias is one, whereas a.t n = 0.6 the predicted mean is nearly 12 times the COBE 
result even for the natural inflation scenario, seemingly requiring a bias of over three. 
However, that conclusion neglects the statistical nature of the quadrupole prediction. 

As discussed above, the quadrupole prediction is not for a unique value but for a 
xt probability distribution function {pdf) with the given mean. For a given observer, 
the prediction is for a given realisation from that distribution. Let us discuss the pre- 
dictions in terms of a quantity q = 10^° 6^ Q^. Then the mean value q = 10^° 6^ {Ql)^ 
and the xHq) pdf with that mean is 

This distribution has a very broad spread (the variance being 2/5 of the mean 
squared), and thus realisations of it may differ considerably from the mean. The 
spread is indicated in figure 5 by the vertical bars through the mean predictions. 

Let us for the time being assume the COBE measurement to be perfect {ie not 
subject to observational errors). Then one can only exclude a given theory on the 
basis of some exclusion level, where the experimental result is far along the tail of 
the distribution. A conventional choice would be that 95% of the distribution should 
predict values above (or below) the experimental measurement. Thus one can exclude 
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at 95% confidence only those theories which predict a mean sufficiently high that 95% 
of the distribution is above the COBE value of g'exp = 0.25 6^. This one can readily 
calculate via error functions, to discover that one can only say with 95% confidence 



that q < 1.16^, corresponding to y {Ql) < 1-05 x 10~^fo. 

One can loosen this constraint even further by incorporating the fact that the 
COBE results possess experimental errors. One can model the COBE data with a 
probability distribution function describing the expected values were the experiment 
to be repeated — a sensible choice might be to assume Q^2^ to be gaussian distributed 
with mean 5 x 10~^ and width 1.5 x 10~^ — and use this to construct a pdf Pexp(?exp) for 
the experimental measurement gcxp- One then tests each member of the experimental 
distribution for exclusion and takes a weighted mean, thus constructing a rejection 
functional on the predicted theoretical pd^ as 



roc rgcxp 
7^ = Pexp(gexp) 1 - 2 XbWq 



dqe^p (4.12) 



where | signals the modulus. Defined on one side of the probability distribution 
like this, a value TZ > 0.9 signals a 95% rejection of the theoretical pdf in the light 
of the modelled experimental data. With the experimental modelling as suggested 
above, then predicted mean values of q up to 2.36^ are allowed, corresponding to 95% 
confidence that y/iOfi < 1.5 X lO-'b. 

Hence one sees that the statistical uncertainties in the theoretical quadrupole 
prediction make it of little use in constraining theories. For example, with this last 
constraint n = 0.6 is allowed with a mild bias of 1.15. 



4.2 The variance at 10" 

Following the filter function formalism of Bond and Efstathiou (1987) (see also Efs- 
tathiou 1991), the mean (over all observer points) of the anisotropy seen by a given 
experiment is given by 

{if))-ii:(^l + mlF, (4.13) 

where Fi is a filter function appropriate to a given experimental configuration. 

For ground-based experiments which typically feature two and three beam con- 
figurations, this filter function can be rather complex. For the COBE experiment 
it is much simpler, as COBE reconstruct the smoothed fiuctuation field across the 
entire sky and calculate the variance of that. The original COBE beam is well ap- 
proximated by a gaussian with Full Width Half Maximum of 7°. However, before 
calculating the variance they smooth again by convolving with a further 7*^ FWHM 
gaussian, a procedure equivalent to an original smoothing by a 10*^ FWHM gaussian. 
Such a gaussian has a variance = (4.25°)^. Thus the appropriate form of the filter 
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function is 



Fi = - exp (- (4.257r//180)') 



(4.14) 



It is now trivial to calculate the predicted 10° variance as a function of n, and 



multiplied by the 



in figure 6 we plot its square root, henceforth denoted ^ 
bias, for both power-law and natural inflation. In the latter case, one sees that the 
results are essentially exactly linear. We have been unable to show why this should 
be analytically. The appropriate fitting function is 
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and so for power-law inflation one has 



(4.15) 
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One can also readily calculate the cosmic variance, giving the spread in values 
about these means that would be measured by differently positioned observers. The 
variance of the is 2{Qf)'^/{2l + 1). For the 10° result, this is 10% at n = 1, rising 
to 12% at n=0.6 (this result remains true whether or not there is a gravitational 
wave contribution to the anisotropics). This is a negligible correction to the larger 
COBE error in the present observations, but is ultimately a hmiting obstacle to the 
conclusions one can draw on this large angular scale. 



4.3 Microwave background constraints 

In estabhshing constraints on the bias at fixed n, it is the 10° data which are of 
primary interest. Some particular values worthy of note are that for the flat spectrum 
- 1.05 X lO^V^, and that for n = 0.6 of ^ = 5.17 x lO'V^ for power- 

IQO ' ^ 100 ' 



off: 



law inflation and ^ „ = 2.99xl0 ^/6for natural inflation. These of course are to 



= ('1.1±0.2) X 10-^ 



be compared with the COBE observations of f- 

The mean quadrupole prediction is also very similar to the flat spectrum predic- 
tion with bias one, but the statistical uncertainties make this comparison rather less 
relevant. It is worth recalling that an extrapolation of the quadrupole from the COBE 
data at smaller angles, assuming a power-law spectrum, gives a somewhat larger pre- 
diction for the mean quadrupole of 16±4/xi^, corresponding to (6±1.5) x 10~^ (Smoot 
et al 1992). Note that this extrapolation estimates the mean quadrupole, not a spe- 
cific realisation. Hence this again favours values of the bias not much exceeding 1, 
unless one inserts large scale power to boost the theoretical prediction for the mean. 

Lacking full access to the COBE data, one must make an operational choice as 
to what to take as the observational limits. The COBE team provide fits of their 
data to power-law spectra, where both the slope and amplitude are treated as fitting 
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parameters. The two pertinent pieces of information are that the observed (root of 
the) variance at 10° is (1.1 ± 0.2) x 10~^, and the allowed slopes in the fit carried out 
by the COBE team are n = 1.1 ± 0.5, where both errors are 1-sigma. These are of 
course not independent, as 10° is at the lower end of the fitting range. Nevertheless, if 
one takes the upper 2-sigma value of the extrapolated quadrupole (24/iii'), then this 
actually gives a prediction for the 10° variance which is well above the 2-sigma limit 
on the 10° data (for n = 1, the root of the mean 10° variance is about 2.03 times the 
mean quadrupole, the relative factor falling as n decreases). This conclusion (also 
implicit in Efstathiou, Bond and White (1992)), argues that the fits are not rigid 
enough to be of much use in constraining theories. Because of the cosmic variance 
and also experimental uncertainties, the quadrupole itself also appears of little use in 
this type of analysis. Consequently, we choose to adopt simply the 10° result as it 
stands, without further incorporation of the COBE data. We relax their error bars to 
2-sigma, and for many purposes we are only interested in the upper limit thus given 
of 1.5 X 10~^. This is a particularly useful way of utilising the results, as it seems 
likely that this number can only go downwards if it is to avoid conflict with other 
experiments and so any limits quoted on its basis are likely only to become stronger 
with improved observations. 

The full limits on the bias as a function of n from this criterion, for both styles of 
inflationary model, are plotted in figure 11 at the end of section |^ along with other 
constraints derived in that section. Some sample results are that for n = 1 we have 
the obvious 0.7 < b < 1.5, while for n = 0.6 we require 3.4 < b < 7.4 (power-law 
inflation) and 2.0 < b < 4.3 (natural inflation). Let us finally recall the uncertainties 
of normalisation. Were one to have used the J3 spectral normalisation rather than 
(Js, these bias values would be lowered by about 10%; the use of the Bardeen et al 
(1986) transfer function rather than the one we use would reduce the required bias 
by a further 10%, and a higher baryonic content raise it by a similar amount. 

5 Intermediate Scales: Galaxy Clustering and Bulk 
Velocity Flows 

Now we consider scales which are big enough that the filtered density contrast is 
still evolving linearly, yet small enough that there is information about the galaxy 
correlation function and about the peculiar velocity field. Thus we are discussing 
scales from around 10/i~^ Mpc up to perhaps 100/i~^ Mpc. The galaxy correlation 
function is sensitive to the combination M(M, x), whereas the peculiar velocity field 
is is sensitive to (5(M, x) itself. 
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5.1 The galaxy correlation function 

There is ever increasing evidence that the amount of large-scale galaxy clustering 
is greater than can be accommodated in the standard CDM cosmogony. The most 
striking piece of evidence is provided by the measurement of the galaxy angular 
correlation function w{9) in the APM survey (Maddox et al 1990, 1991), based on a 
sample of over two million galaxies. Further evidence pointing to the same conclusion 
has been provided by the 'counts in cells' of the QDOT survey (Saunders et al 1991), 
and more recently by a survey of redshifts of APM galaxies (Dalton et al 1992), and 
in the power spectrum inferred from the CfA survey (Vogeley et al 1992) and that 
inferred from the Southern Sky Redshift Survey (Park, Gott & da Costa 1992). 

Our earlier work with Will Sutherland (LLS) on structure from the power-law 
inflation model was primarily concerned with an investigation of whether or not this 
excess clustering could be intrinsic to the primeval spectrum, were the spectrum to 
be of power-law form. In this subsection we recap briefly on that analysis and its 
conclusions. The remaining results of that paper, concerning microwave anisotropics 
and the formation of nonlinear structure, are superceded by the other sections of this 
paper. 

On the large angular scales where the discrepancy arises, the angular correlation 
function w{9) can readily be calculated in linear theory. It is given (Peebles 1980) by 
Limber's equation as an integral over the galaxy correlation function C,gg{r) = &^^(r), 
where the mass correlation function ^{r) is defined in section ^. A linear calculation 
is expected to be reliable on scales above about 2°. The results appear in figure 
7, reproduced from LLS. The fiat spectrum n = 1 falls well below the observational 
data, but one can see that as the spectrum is tilted, the extra power does indeed make 
itself evident in the clustering statistics. It is suspected that the APM survey may 
contain small residual systematics which bias the observational estimates upwards, 
so it seems reasonable to regard values of n from 0.3 to 0.6 as good fits to the excess 
clustering data. (This is actually rather more restrictive than the range allowed by 
Efstathiou, Bond and White (1992), as we see below.) We emphasise once more that 
this result is independent of the value of the bias parameter. 

It has become a common reference point to consider spectra with low effective flh 
when considering clustering data. Figure 8 compares the processed spectra for n = 1, 
0.6 and 0.5 with that of a CDM model with a fiat spectrum but with Qh ~ 0.2. Such a 
model has been much touted, especially by Efstathiou and collaborators (Efstathiou, 
Sutherland & Maddox 1991; Efstathiou 1991; Efstathiou, Bond & White 1992), as 
providing a good fit to all of the excess clustering data. One sees that the tilted 
power-law spectra resemble this spectrum much more closely than does standard 
CDM, though with some additional reduction in short scale power. The Qh ~ 0.2 
model has also been claimed as a good fit to the distribution of rich clusters (Lilje 
1992; Scaramella 1992), to the CfA survey (Vogeley et al 1992) and to the Southern 
Sky redshift survey (Park, Gott & da Costa 1992), all of which appear inconsistent 
with standard CDM. We thus expect tilted spectra to also do well on these criteria. 
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though we have not attempted a direct comparison here as we lack suitable data. 

One can quantify this connection somewhat by utilising a quantity introduced 
by Wright et al (1992) called the excess power, which is a functional of the power 
spectrum defined as 

r^n a(25/i"iMpc) , , 

EVP] = 3.4 \ , 5.1 

For standard CDM and our transfer function, E = 0.95 (there is a very weak depen- 
dence on the choice of transfer function, which still plays a role on these scales). We 
find a reasonable fit is given simply by 

E[n] = lAA-nj2 ; 0.3 < n < 1 (5.2) 

This can be related to the F parameter of Efstathiou, Bond and White (1992) as 

r = if^V°'' (5.3) 



2 V2.88-n 

The range of F these authors consider a reasonable fit to the APM data is 0.15 < 
F < 0.30, corresponding to 0.15 < n < 0.67, which we see is actually rather looser 
than the range we took above. 

The relation between n and F provides an excellent means to compare results of 
different authors on intermediate length scales. It is worth remembering though that 
well away from the scales 8/i~^ and 25h~^ Mpc used to define excess power the spectra 
do show significant differences, and so this relation is not as useful for considering 
the microwave background anisotropies, and even less so on short scales where the 
power-law spectra show considerably less power than their F-equivalents. 

It would be invidious of us to close this subsection without reference to a huge 
body of literature which explains the excess clustering by a variety of means other 
than by adding intrinsic power to the primeval spectrum. One way in which this 
can be done is to modify the transfer function. Introducing a cosmological constant 
with Q\ ~ 0.8 provides a realisation of the F ~ 0.2 standard CDM model discussed 
above (Efstathiou, Sutherland & Maddox 1990; Efstathiou, Bond & White 1992), 
which remains viable though under threat from observations of quasar lensing (see eg 
Fukugita et al 1992). A 17 keV neutrino with a lifetime of the order of 10 years can 
also provide an appropriate modification of the transfer function (Bond & Efstathiou 
1991), though the evidence for such an object from particle physics seems to be 
evaporating. Hybrid models combining both hot and cold dark matter in carefully 
chosen combination are also finding increasing favour (Schaefer & Shafi 1992; Davis, 
Summers & Schlegel 1992; Taylor & Rowan- Robinson 1992). 

An alternative strategy is to modify the galaxy formation process itself. A partic- 
ularly elegant proposal is to allow quasars to suppress local galaxy formation (Babul 
& White 1991), though this must now contend also with the COBE measurements 
implying a bias of around 1. Cooperative galaxy formation, a phenomenological (at 
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present) model in which galaxy formation is favoured in the neighbourhood of other 
galaxies, has also been proposed (Bower et al 1992), which can have the effect of a 
scale-dependent bias. Finally, it has been suggested that a detailed analysis of nonlin- 
ear evolution assuming a bias of around 1 (made more attractive by COBE) coupled 
with a velocity antibias shows that there is no conflict between the predictions of this 
model and the observed galaxy clustering (Couchman & Carlberg 1992). 

5.2 Bulk peculiar velocities 

There are many controversial aspects to the measurements of bulk velocity flows, but 
nonetheless they provide a useful measure of the absolute magnitude of the power 
spectrum on intermediate scales (for a recent review, see Kashlinsky & Jones (1991)). 
Quite detailed work on the bulk velocity flows has been carried out for the power- 
law spectrum (and also allowing Q < 1) by Tormen, Lucchin and Matarrese (1992) 
and by Tormen et al (1992), though our analysis will take a somewhat different and 
simplified approach. 

Just as one can filter the density contrast, one can also filter the peculiar velocity 
field to obtain a bulk flow quantity, on any desired scale. As long the density contrast 
on that scale is in the linear regime, the bulk flow can be constructed from the linearly 
evolved Fourier components Vk, given by (Peebles 1980; Kolb & Turner 1990) 

Vk = ^k— 5k (5.4) 

where k is the unit vector in the k-direction. 

Through the galaxy redshifts, the peculiar velocity field at the present epoch is 
observable within a sphere around us whose radius is a few hundred Mpc. Smearing 
over a scale > Wh-^ Mpc, one obtains a bulk velocity field which is still evolving 
linearly. Taking the mass density contrast as being times the galaxy number 
density it too can be observed, and inserting both quantities into Eq. (|5.4| ) the value 
of b can be estimated. This method of estimating b has been extensively explored for 
the standard case n = 1. Until recently, applied for instance to the region around the 
Virgo supercluster, it seemed definitely to require a bias factor b significantly bigger 
than 1 (Peebles 1980; Kolb & Turner 1990). With the advent of better measurements 
this is no longer the case, and a value 6 = 1 is apparently possible (Dekel 1991). 
Without direct access to the data it seems difficult to know how this result is affected 
if n < 1. The possibility of a velocity bias (Carlberg, Couchman & Thomas 1990; 
Couchman & Carlberg 1992) further complicates the analysis. 

We base our comments here around the recent results concerning IRAS galaxies 
form the QDOT survey. Our analysis essentially mimics that of Efstathiou, Bond 
and White (1992). One gets two useful statistics from the QDOT survey. The first is 
that, in combination with bulk flow data, one can estimate directly the bias of IRAS 
galaxies. In general, the bias of the infra-red selected (and thus typically young) 
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IRAS galaxies will not be the same as that of the optically selected galaxies discussed 
thus far, and so we denote this bias by bj. It is often stated that IRAS galaxies 
are somewhat less clustered than their optical counterparts; for example Saunders, 
Rowan-Robinson and Lawrence (1992) suggest bj = (0.69 ± 0.09)6 at the 1-sigma 
level. According to Taylor and Rowan- Robinson (1992), there are three independent 
dynamical estimates of 6/, giving respectively bj = 1.23±0.23, 1.16±0.21 and 1.2±0.1 
where all errors are 1-sigma. One cannot statistically combine errors from the same 
data set, so we adopt as a result (which we assume has something like 2-sigma errors) 
that bi = 1.2 ±0.3. 

The second useful statistic is the dispersion of counts in cells in 30h~^ Mpc cubes, 
given as 0.46 ± 0.07 (Saunders, Rowan- Robinson & Lawrence 1992) where we have 
doubled the quoted error bars. Because we have a good estimate of the bias, we can 
immediately calculate the mass variance in these cubes, and then use the spectrum 
to calculate the variance at any other scale. In particular, one can scale down to 8h~^ 
Mpc and thus determine limits on the bias parameter b. Note that the use of the 
optical bias parameter is just a way of representing the results, without implying any 
reference either to optically selected galaxies or to scales on the point on nonlinearity. 
We note also that even if one did not know the variance of the QDOT counts in cells, 
one would already have strong limits from the knowledge that optical galaxies are 
not too much more clustered than IRAS ones. 

The comparison between the bias limits obtained by this method and those ob- 
tained from the COBE measurements represent a particularly strong combination. 
In our earlier paper (Liddle & Lyth 1992) we pointed out that because of the strong 
gravitational contribution to the microwave anisotropics, this comparison rules out 
at least the simpler models of extended inflation because of the n < 0.75 bubbles 
constraint on these models. In figure 11 at the end of this section, these constraints 
appear together and highlight the allowed regions of the n-b parameter plane. 

Before finishing this section, we mention one other simple comparison which can 
be made — to the bulk flows in spheres as measured by POTENT. This approach 
considers the dispersion of = |v|, given on each scale through the spectrum of v, 



Using this approach it is easy to consider the effect of n. Figure 9 illustrates the scaling 
of the dispersion with i?j for a top hat fllter. Smaller values of n give signiflcantly 
higher predictions, but the scaling with Rf is roughly the same. 

To compare this result with observation we use the velocity fleld reconstruction 
results from the POTENT method, pioneered by Bertschinger, Dekel and collabora- 
tors (Bertschinger & Dekel 1989; Dekel, Bertschinger & Faber 1990; Bertschinger et 
al 1990; Dekel 1991). It is weU known (Kolb & Turner 1990) that typical theories, 
including standard hot and biassed cold dark matter models and also models seeded 
by topological defects, tend to predict bulk velocities rather lower than those ob- 
served, particularly for standard CDM with high bias. Nevertheless, it is not trivial 




(5.5) 
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to compare theory with observations, because one predicts only the rms velocity aver- 
aged over the entire universe and the observed distribution may not be a fair sample. 
In particular, the possible contaminating effects of the great attractor (if backside 
infall can be unambiguously identified (Mathewson, Ford & Buchhorn 1992)), may 
distort observations. In the Monte Carlo work analysing bulk flows of Tormen et al 
(1992), special criteria for choosing the observer points are employed before statistical 
comparisons are made. 

The POTENT results are obtained via a two-stage smoothing (Dekel 1991). First, 
the original data is smoothed with a gaussian of radius 12h~^ Mpc, and then a top 
hat of radius -R/ is used to provide the published data. In figure 10, we apply this 
two-stage smoothing to our spectra, which reduces the predictions, especially at short 
scales. These are compared with the POTENT data at different biasses. It is seen 
that standard CDM at bias 1 performs rather well, significantly better than n = 0.6 at 
bias 2. Nonetheless, it appears that either model is still tenable given the statistical 
nature of the quantity being observed. In particular, the influence of large scales 
filters down significantly to small scales with this statistic, so an unexpectedly large 
fluctuation (and the great attractor would certainly be one in these CDM models) 
can easily move the whole observed data set to well above the mean, as predicted 
from the spectrum, that would be observed in a fair sample. 

Some numerical analysis (based on Monte Carlo simulations) of bulk velocity 
flows with power-law spectra has been made by Tormen et al (1992), who consider 18 
models combining all possible combinations of 6 = 1, 1.5, 2; n = 1, 0.5, 0; = 0.4, 
1. Of most interest to us is the Q = 1, n = 0.5, b = 1.5 model, which most closely 
represents the criteria we are establishing in this paper. They employ a maximum 
likelyhood test, and find that this model is twice as likely as the model with flat 
spectrum and bias 1.5, and is four times more likely than the flat spectrum with bias 
1. 

5.3 Summary of the intermediate/large scale results 

One of the key aims of this paper is to investigate the constraints in the n-b parameter 
plane. The most significant constraints from large and intermediate scales are plotted 
in figure 11, and all the lines on this figure are to be taken as representing the 2-sigma 
range on either side of the (unplotted) mean value. The constraints plotted are 

• The upper limit on n (independent of b) obtained from requiring a fit to the 
APM data (dotted line). We have taken our constraint to be < 0.6 (see figure 
7), though we note that Efstathiou, Bond and White (1992) allow a range for 
their F which is roughly equivalent to letting n be as large as 0.67, which is very 
conservative. 

• The limits on bias as a function of n from the QDOT survey (dot-dashed lines), 
as discussed in section ^2. In a sense these are the weak link in the diagram. 



23 



as although we have been very conservative in the errors on this observation it 
is not inconceivable that changes will be seen in the future. Nonetheless, we 
believe that these limits should be taken very seriously and may well strengthen. 

• The COBE limits on bias as a function of n for power-law and extended infla- 
tion (solid lines) and for natural inflation (dashed line), as discussed in section 
^.3. These indicate the upper and lower 2-sigma range of the COBE 10° data. 
These are very strong results, as it seems inevitable that if any change to the 
observations occurs it will be to lower the fluctuations and thus increase the 
required bias. 

By themselves, the QDOT and COBE data allow a region of the n-h plane which is 
roughly a triangle with corners {n,h} = {.7, 1.6}, {1, 1.6} and {1, 1}. Leaving aside 
for the moment the problematical galaxy clustering represented by the APM data, 
we ask in the next section whether data on smaller scales can discriminate between 
points within this triangle. 

6 Small Scales: Galaxies and Galaxy Clusters 

Now we study small scales < lO/i^^ Mpc, corresponding to M < IO^^Mq. This 
is the range on which gravitationally bound objects exist, in particular galaxies with 
mass IO^Mq < M < IO^^Mq, and galaxy clusters (including the small ones known as 
groups) with mass lO^^M© <M< IO^^Mq. 

In contrast with the treatment so far, we shall not attempt in this section to 
work at the 2-sigma level. The uncertainties in observational data, where they are 
considered at all, will be taken simply as those quoted in the literature. 

For small scales, there are many quantities which one might hope to calculate and 
compare with observation. A partial list is the following. 

• The redshift of formation of a given class of objects, such as bright galaxies or 
rich galaxy clusters. 

• The number density n{z, > M) of all gravitationally bound systems with mass 
bigger than M, which exist at redshift z. 

• The mean virial velocity of a given class of objects. 

• The number density n{z, > v) of all gravitationally bound systems with virial 
velocity bigger than v, which exist at redshift z. 

• The correlation function 'Cobj('") of a given class of objects. If it has the same 
shape as the correlation function ^(r) of the mass over scales of interest, it is 
useful to define a bias parameter by ^obj = &obj^- 
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• The angular correlation w{6) of bright galaxies, on the angular scales 6 < 2^ 
which correspond to linear scales r < 10h~^ Mpc. 

• The dispersion a^{r) of the relative velocity of a pair of galaxies separated by 
distance r. 

As galaxies and clusters are gravitationally bound objects, these quantities can- 
not be calculated using only the linear cosmological perturbation theory which was 
adequate for intermediate and large scales. Considerable insight can, however, still be 
gained from the linear theory if it is combined with order of magnitude estimates in- 
spired by a spherical model of gravitational collapse. The general conclusion from this 
approach is that the n = 1 case is viable, with b somewhere in the range 1 < 6 < 2. 
In the context of the present paper, we wish to ask whether the linear approach will 
also permit other values of n and b, within the triangle allowed by the QDOT and 
COBE data. 

In the n = 1 case, small scales have also been studied through non-linear calcula- 
tions, namely A^-body simulations. Most such calculations agree that the n = 1 case 
is viable on small scales (Davis et al 1985; Gott et al 1986; White et al 1987; Carlberg 
& Couchman 1989; Frenk et al 1990; Carlberg et al 1990; Bertschinger & Gelb 1991; 
Couchman & Carlberg 1992), though some find that it has problems (Suto et al 1992, 
and other references sited there). A useful summary of the situation for n = 1 has 
been given by Davis et al (1992). The corresponding calculations for n < 1 are so 
far in their infancy (Vogeley et al 1992; Park et al 1992; Cen et al 1992), but so far 
they seem to suggest that a value n < 1 is preferred. We will not consider numerical 
simulations further in the present paper. 



6.1 Gravitational collapse 

The linear approach makes essential use of the filtered density contrast 5(M, x). Fol- 
lowing the usual practice for theoretical calculations, we use the Gaussian filter, de- 
fined by Eqs. (|3.11| )- (|3.14|) . the dispersion of the filtered density contrast is thus 



a(M)= V{k)exp{-R}e)^ (6.1) 

Jo fb 

At this point we need to distinguish between the density contrast of the cold dark 
matter, and that of the baryonic matter. They are the same on filtering scales M 
well in excess of 10^ Mq. As M is reduced below this figure, the density contrast 
of the baryonic matter falls off sharply, because pressure forces prevent its growth 
until long after decoupling (Peebles 1980; Kolb & Turner 1990). In contrast the cold 
dark matter has negligible pressure, and its density contrast grows like (1 + z)~^ just 
as it does on larger scales. This difference is presumed to account for the absence 
of luminous galaxies with mass below IO^Mq. In what follows, 5(M, x) will always 
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denote the density contrast of the cold dark matter, though we are usually interested 
only in scales bigger than IO^Mq. 

As long as it is less than 1, the filtered density contrast evolves linearly except in 
those rare regions of space where it exceeds 1 in magnitude. In the 50% of these regions 
where it is positive, gravitational collapse takes place. At least initially, the collapsing 
regions have mass > M, because the filtered density contrast does not 'see' structure 
on smaller scales. [] In the approach that we are considering, gravitational collapse is 
modelled by taking the collapsing region be spherically symmetric. One then finds 
(Peebles 1980, Eqs. 19.50 and 19.53) that when a given mass shell stops expanding, 
the mean density inside it is a factor 97r^/16 times bigger than the mean density of the 
universe. If the matter inside the shell had evolved linearly (density contrast oc (1 + 
z)~^), its density contrast at that time would have been S = (3/5)(37r/4)^/^ = 1.06. 
After it has stopped expanding, the shell collapses. If spherical symmetry continues 
to hold, and one neglects pressure forces, the shell has collapsed to a point by the time 
that the age of the universe has increased by precisely a factor 2, corresponding to a 
mean density contrast within the linearly evolved shell of (5 = (3/5)(37r/2)^/^ = 1.69. 
Numerical studies (eg. Peebles 1970) indicate that by about this time, pressure 
forces will in fact have virialised the random motion of the constituents of the object. 
After this initial virialisation, the object can lose energy (dissipation) which further 
increases its virial velocity. 

The use of the spherical collapse model in the present context is somewhat prob- 
lematical. First, filtering on a mass scale M will distort the profile of a peak with 
mass of order M, making it lower and broader. Thus the linearly evolved filtered 
density contrast will not have precisely the behaviour described above, even for a 
spherical peak. Second, the departure from sphericity of a linearly evolved collapsing 
region 'seen' by the filtered density contrast can be estimated from the linear theory. 
As we shall see it is quite significant, and that of the true (unfiltered) density contrast 
will be more so. Finally, any initial departure from sphericity is amplified during the 
collapse (Peebles 1980). Nevertheless, the fact that the peaks of the density contrast 
are at least roughly spherical (during the linear regime), as well as numerical sim- 
ulations that have been done for the n = 1 case, indicate that two features of the 
spherical collapse model translate into reality. 

• The regions with mass > M which have undergone gravitational collapse can 
be at least approximately identified with the regions where the linearly evolved 
density contrast 5(M, x) exceeds some threshold 6c- 

• A collapsing region does not fragment into a large number of separate objects, 
which means that the mass of the resulting gravitationally bound systems is 
also > M. 

The optimal choice of the threshold 6c is a matter of debate. Many authors take 
^As discussed later, this statement needs modifying when the collapsing regions are very rare. 
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the value 1.69 inspired by the above spherical collapse model. On the other hand, 
comparison of the hnear estimate of n{> M) (described below) with the estimate 
from numerical simulations suggests a smaller value, Carlberg and Couchman (1989) 
advocating 5c = 1.44 and Efstathiou and Rees (1988) advocating Sc = 1.33 (but 
see also Brainerd and Villumsen (1992)). When making numerical estimates we set 
6c = 1.33, while keeping in mind the effect of taking a bigger value. 

As the collapsed regions 5{M,x) > 5c represent exceptionally large fluctuations 
of a Gaussian random field, there are powerful mathematical results concerning their 
stochastic properties. We shall use some of them in what follows, drawing extensively 
on the work of Bardeen et al (1986, henceforth BBKS). 

6.2 The epoch of structure formation 

On each mass scale, the linear regime ends at the epoch z^\{M) when the filtered 
linearly evolved dispersion a{z,M) of the density contrast is equal to 1. Since this 
quantity is proportional to {1 + z)~^, the epoch is given by 

l + z„i(M) = (7o(M) (6.2) 

where the subscript denotes the present value. At this epoch, the formerly rare 
regions where gravitational collapse is taking place on scales ^ M become common. 
Before studying the linear regime which is our main concern, let us ask how this 
gravitational collapse process proceeds. In particular we want to know whether the 
bottom-up picture of structure formation, known to occur if n = 1, will also occur if 
n<l. 

At the epoch z-ai{M) the collapsing regions, which have mass ^ M, give rise to 
gravitationally bound structures with mass ^ M provided that they do not fragment 
into a large number of pieces. If a{M) is increasing significantly as the mass decreases, 
then at this same epoch the density contrast filtered on mass scales substantially 
bigger than M is still evolving linearly in most parts of the universe. In that case, 
only a small fraction of the mass of the universe is bound into objects with mass much 
bigger than M . Moreover, iV-body simulations indicate that fragmentation does not 
then occur, so that the objects formed at the epoch z^\{M) do not have mass much 
less than M either. The result is a bottom-up scenario, in which gravitationally 
bound systems of successively bigger mass M form at the successively later epochs 
z-a\{M). After a system has formed, it may become part of a bigger gravitationally 
bound system, remaining a discrete object, but it may also merge with other systems. 

The scale dependence of a{M) depends on that of the spectrum V{k), on the 
corresponding scale ~ Rf. In the cosmologically interesting range M > IO^Mq, 
these scale dependences were already shown in figures 3 and 4 both for the standard 
choice n — 1 and for some smaller values. One sees that for n = 1, a{M) increases 
steadily as M decreases, corresponding to the fact that V{k) does not turn over 
as decreases. In contrast, for n < 1 the spectrum V{k) turns over below some 
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mass scale, so that a{M) levels out. This difference in the small scale behaviour arises 
from the small scale behaviour V{k) oc T'^{k)k^^^^\ According to the parametrisation 
Eq. (|3.4| ), T'^{k) is proportional to A;~^ which makes V{k) oc k^~^. This means that 
(t(M) increases logarithmically as M decreases if n = 1, whereas if n < 1 it approaches 
the unfiltered value a, the difference going to zero like (Strictly speaking the 

above statements are modified by logarithmic factors because goes like log^ k 
(BBKS), but they do not change the qualitative behaviour.) 

The small scale behaviour described in the last paragraph is expected to persist 
down to some very small scale, below which the spectrum cuts off sharply and a{M) 
becomes practically equal to the unfiltered quantity a. This scale is called the coher- 
ence scale and depends on the nature of the cold dark matter. If the cold dark matter 
consists of subnuclear particles, the coherence scale is determined by their interac- 
tions. For instance, if it consists of the lightest supersymmetric particle the coherence 
scale is the Hubble scale at the epoch when the particle becomes non-relativistic. If 
it consists of the axion, the coherence scale is the Hubble scale at the epoch when the 
axion acquires mass through QCD effects. 

On the basis of these considerations, we arrive at the following picture of structure 
formation, which covers the case n < 1 as well as the familiar case n = 1. The first 
structure forms at some epoch given by 

1 + Znl = CTq (6.3) 

where ctq is the unfiltered quantity. For n = 1, this epoch is very early and depends 
on the nature of the cold dark matter. Decreasing n makes the epoch later, and to 
an increasingly good approximation allows it to be calculated without knowing the 
nature of the cold dark matter. After this epoch, structure forms according to the 
bottom-up picture. 

Let us estimate the maximum mass Mmax, above which structures collapsing at 
a given epoch Zn\{M) are rare. It should be such that o"o(Mmax) is significantly less 
than 1, so that the density contrast filtered on the scale Mmax is still evolving linearly 
almost everywhere in the universe. To be definite, let us require that it is still evolving 
linearly in 90% of the volume of the universe. From the Gaussian distribution, this 
corresponds to c'"o(Minax) = 1/1.63 = .61. For M = lO^M© this criterion gives 
Mmax = IO^Mq (IO^^Mq) for n = 1 (.7). For this value of M, the epoch 2ni(M) is 
a slowly varying function of M, so the conclusion is that a broad band of masses 
downwards of M = lO^M© to IO^^Mq collapse at the same time, around the epoch 
z^i{10^Mq). For M = lO^^M©, the criterion gives M^ax = 2 x lO^^M©, more or less 
independent of n. For this value of M, the epoch Zni(M) is a fairly rapidly varying 
function, so the conclusion is that a fairly narrow band of masses M = lO^^M© to 
about IO^^Mq collapse around the epoch Zn^lO^^M©). 

We conclude that the bottom-up picture is more or less the same for any n in the 
range .7 < n < 1. A broad range of masses below about lO^M© collapses at about 
the same time, but subsequent collapse takes place in an increasingly narrow mass 
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range. 

The epoch Zni{M) can be read off figure 4, for n = 1 and n — .7. It becomes more 
recent as M or n are reduced. Let us consider three representative values of M. First 
take M = lO^M©. Then ^ni(M) is quite early, 1 + ^ni(M) = 18/6 (9/6) for n = 1 
(.7). A broad range of mass scales downwards of about IO^Mq (IO^^Mq) collapse at 
around this epoch. Second, take M = 10^^ Mq. Then 1 + Zni{M) = 4.5/6 (3.7/6) for 
n — 1 (.7). This is the epoch when a significant fraction of the mass of the universe 
collapses into objects with mass of order IO^^Mq, the mass of large galaxies. However, 
the favoured explanation of a bias factor 6 > 1 is that luminous galaxies originate 
from exceptionally high peaks of the evolved density contrast, and therefore form well 
before the epoch Zni{M) (BBKS; Efstathiou 1990). Finally consider M = IO^^Mq. 
Then, l+2;ni(M) = .82/6, almost independent of n. This means that the epoch Zni{M), 
at which a significant fraction of the mass collapses into large galaxy clusters, hes in 
the future. The prediction is therefore that the observed large clusters originated from 
exceptionally high peaks of the density contrast, which again implies a bias factor for 
these objects (Kaiser 1984). The bias factor for galaxies and galaxy clusters will be 
discussed quantitatively below. 

Finally, we consider a different criterion which is sometimes proposed for the epoch 
of structure formation, namely that it is the epoch when the linearly evolved spectrum 
V{k) is equal to 1 on the corresponding scale k^^ = Rf. This is practically equivalent 
to the criterion a{M) = 1 if V{k) is significantly increasing, say like /c™ with m ~ 1, 
because then the filtered spectrum is 'P{k) exp(— /c^i?^) ~ V{Rj^)6{log{kRf)), where 
S denotes the Dirac delta function. The condition for this to be the case is more or 
less the same as the condition for bottom-up structure formation to occur. It fails, 
in particular, when k corresponds to the maximum of V which occurs when n < 1, 
so it is not true that in that case the first structure forms at the epoch when 1 + 2; is 
equal to the maximum value of V{k), as was incorrectly stated in LLS. 



6.3 The statistics of the collapsed regions 

Now we study the stochastic properties of the collapsed regions, defined as regions 
where the hnearly evolved density contrast 5{M, x) exceeds a threshold Sc- In them, 
5(M, x) is more than u standard deviations above zero, where 

iy{z,M) = 6Ja{z,M) (6.4) 

= Sc{l + z)/ao{M) (6.5) 
. 1 + z 



^l + ^nl(M) 



(6.6) 



We are working in the linear regime, corresponding to > 6c > 1- 

The collapsed regions occupy a volume fraction V given by the Gaussian distri- 
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e-"/-" (6.7) 



but ion, 
leading to 

V{u) = erfc(z//V2)/2 (6.8) 
= {27T)-^/^u-^e-''/\l - u-^ + 0{u-^)) (6.9) 

For u = 1, 2, 3, 4 the volume fraction is V = .16, .023, .0013, .000031. In practice one 
is not interested in values > 4, because the collapsed regions are then too rare to 
be physically significant. The corresponding mass fraction is about (1 + 6c) = 2 to 3 
times bigger than the volume fraction. 

To say more one needs to know the shape of the spectrum. We shall list the 
relevant results given by BBKS. They involve only two moments of the spectrum, 
defined by 



< k\M) > = a-\M) / k^exp{-eR})V{k)-^ (6.10) 

< kHM) > = a-\M) / k''eM-k^R))Vik)^ (6.11) 

Jo ■' k 



The quantity (/c^) is the mean of the operator, ie., of the quantity 5~^V^(5. Simi- 
larly, (/c^) is the mean of V^. 

A relevant length scale is defined by Rl = 3{k'^)/{k^). For any spectral index 
n > — 1, it is easy to show that in the limit of small filtering scale Rf, 

R / 6 

* ^ ^ (6.12) 



Rf \l + n 

For the case of CDM with .7 < n < 1, the ratio is in the range 1 to 3 for the entire 

range of cosmologically interesting masses. 

Another relevant length scale is (/c^)^^/^. On large filtering scales, such that V{k) 
is increasing fairly strongly at k~^ ~ Rf, the ratio is close to 1. As 

the scale is reduced it increases, but is <^ 10 in the cosmologically interesting range 
M > IO^Mq. 

Finally, it is convenient to define the dimensionless parameter 

7(M) = {e)/{kY' (6.13) 

It falls from about .7 to about .3 as M decreases from 10^^ Mq to IO^Mq, for .7 < 
n < 1. 

For sufficiently large u, each collapsed region is a sphere surrounding a single peak 
of 5. However, the departure from sphericity is considerable in the cosmologically 
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interesting regime. BBKS show that a quantity x ^, which is roughly the fractional 
departure from sphericity, is well approximated by 

a; = 71/ + 6* (7, 71/) (6.14) 

where 

^ 3(1-0^^) + (1.216 -.9o^«)exp|---,/2(W2)-] 

|3( 1 - 1-2) + .45 + (7i'/2)2] + 

The asphericity is plotted in figure 12 for M = lO^M© and 10^"^ Mq, for both n = 1 
and n = .7, and is seen to be ^ .3 even at i/ = 4 and IO^^Mq. We emphasise that 
this is the asphericity seen in the linearly evolved, filtered density contrast. The 
asphericity in the true, unfiltered density contrast will be bigger, and will increase 
during collapse. 

Three useful number densities arc given by BBKS. First, the density of the 
Euler number of the surfaces bounding the collapsed regions is 



Second, the number density of upcrossing points on these surfaces is 
nnp{i^, < A;^ >,7) 



57^(1 -V/9)V- ^^^^-^"^^^/^^^. 



(6.16) 



-1/2/2 



(6.17) 

An upcrossing point on a surface of constant 6 is defined as one where V5 points 
along some chosen direction. 

The third number density is ripeak, the number density of peaks which are more 
than u standard deviations high. BBKS give expressions for npeak, but they point 
out also that in the cosmologically interesting regime it is quite well approximated by 
n-ap. We shall use this approximation in what follows. It suggests that if a collapsed 
region contains several peaks, they are not buried deep inside it; rather, the boundary 
of the region is presumably corrugated, wrapping itself partially around each peak. 

In the limit i/7 ^ 1, ripeak = n^p — I^t-x- This is in accordance with the fact that 
in this regime, each surface is a sphere surrounding a single peak. It contribTitcs +1 
to the number of upcrossing points, and +2 to the Euler number. As i/ decreases 
the surface becomes deformed, but at first its contributions to riup and are not 
affected. Eventually though, it may become so corrugated that it has more than one 
upcrossing point, and may become a torus so that its Euler number is less than 2 
(of course it then has more than one upcrossing point whatever its shape). As a 
result, its contribution to becomes less than its contribution to nup. The ratio 
2npeak/^x (approximated as 2n-ap/n^) is plotted in figure 13. (We are not interested in 
the regime i/ < 1, but for the record is negative there, indicating that the surfaces 
5(M, x) = va definitely do not have spherical topology; in fact they percolate leading 
to a sponge-hke topology for the regions (5(M, x) > 5c (Melott 1990).) 
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One would like to know the number density ncou of the collapsed regions. It 
satisfies the inequality < ricow < n^eak^ and so is equal to in the limit 

z/7 ^ 1. According to figure 13, this defines it within a factor 2 for i/ > 3. The 
average number N of peaks per collapsed region satisfies 1 < < 2npeak/'^x' 
according to figure 13 it is no bigger than 2 for > 3. 

6.4 The mass of the collapsed regions 

As the filtered density contrast does not contain structure on scales much less than 
the filtering scale Rf, we expect the average radius of a peak to be > i?/. It can be 
estimated from the number density of peaks with arbitrary height npeak(~oo), which 
is equal to .016i?~^ where the length R^, was defined at the beginning of the last 
subsection. This is also the number density of troughs, so a rough estimate of the 
average radius of a peak or trough at half height is (2npeak(— oo)"^/^/4 = .8R^. As we 
discussed earlier, R^/Rf ~ 1 to 3 for all cases of interest, so the average peak size is 
roughly of order Rj. Since there are few peaks with radius less than Rf, this suggests 
that the probability distribution of peak sizes is fairly narrow, most peaks having a 
radius around the average. 

We would like to compare the average peak radius Rf with the average radius of a 
collapsed region. Equivalently, we would like to compare the filtering volume Vf with 
the average volume of a collapsed region. The latter is is equal to the volume fraction 
occupied by the collapsed regions, Eq. ( |6.9|) , divided by their number density ricoii- In 
general we do not know ricow, but we do know the peak number density npeak? so we 
can calculated the average volume Vpeak per peak of a collapsed region. It is useful to 
give the result as a fraction to the filtering volume Vf, 



The corresponding mass fraction is roughly Mpeak/M = (1 + 5c) Vpeak/ V/. The volume 
fraction is plotted in figure 14, and one sees that except for large u it is bigger than 
1. 

At large u, the ratio may be calculated using the approximation riup = "^x/^ 
together with Eq. (|6.9|), 



The fact that it is small indicates that a typical collapsed region is sitting on top of 
a peak, as opposed to the situation for smaller u where it encompasses most of the 
peak. This is natural, because for large u a collapsed region is exceptionally high. 




(6.18) 



Vpeak/V; ^ i2nf\{e)/3r'/'u-yVf 



(6.19) 




(6.20) 
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6.5 The number density n(> M) 

The main application of these resuhs is to estimate the number density n{> M) 
of gravitationally bound systems with mass bigger than M, at a given epoch before 
Zni{M). The systems are supposed to be identifiable by looking at the linearly evolved 
density contrast S{M, x). Each collapsed region, defined as one in which S{M, x) > Sc, 
is supposed to contain one or more systems with mass bigger than M. As we have 
just seen this cannot be correct for large u, but let us accept it for the moment. 

If each collapsed region is identified with a single system, then n(> M) = ricoii- In 
general this recipe is useless for lack of an expression for ricoii- A different prescription, 
which does lead to a calculable expression, is to identify each peak within a collapsed 
region with a different collapsed object, 

n{> M) = Upeak ^ Uy^p (6.21) 

This estimate (usually without the simplifying second equality) is widely used in 
the literature. It is certainly the same as the estimate n(> M) = ricoii for large u, 
where we know that there is just one peak per collapsed region. To what extent the 
prescriptions are the same for lower v is not known, because the number of peaks per 
collapsed region is not known. 

If at some epoch the linearly evolved density contrast does have many peaks 
within a collapsed region, an interesting situation arises. At a somewhat earlier epoch, 
5(M, x) was smaller, and a separate contour 5(M, x) = 6c was wrapped around each 
peak. In other words, each peak of the linearly evolved density contrast, filtered 
on scale M, was inside a single collapsed region, and presumably represented a sep- 
arate gravitationally bound system. At the later epoch when the collapsed region 
encompasses many peaks of the linearly evolved density contrast, we have a bigger 
gravitationally bound system. If the original systems survive, the identification of 
each peak with a separate system is correct, but it misses the larger system which 
contains the original systems. Of course, missing this one system does not affect the 
total count much, so if this case is typical of collapsed regions containing many peaks 
the estimate n{> M) = ripcak is better than the estimate n{> M) = ricoii- If, on the 
other hand, the original systems have merged, that identification is wrong, and the 
whole of the collapsed region should be identified with just one gravitationally bound 
system. If this case is typical, the estimate n{> M) — ncou would be better, if only 
we had a formula for it. Which case is the more likely? A clue is provided by the 
observation made earlier, that if there are several peaks in a collapsed region they 
typically seem to lie near the surface of the region, a part of the surface wrapping it- 
self around each peak. This picture would suggest that the estimate n(> M) — ripeak 
is the more reasonable, the peaks of the linearly evolved density contrast in a typical 
collapsed region representing structures which have not existed long enough to merge. 

A different prescription was used by Press and Schechter (1974), to derive a widely 
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used alternative formula. They worked with the differential number density, 



dn _ d 
dM ^dM 



n{> M) 



(6.22) 



At a given epoch, if the filtering mass M is increased by an amount dM then v 
is increased by an amount dv, and the volume fraction occupied by the collapsed 
regions is reduced by an amount dV given by Eq. (|6.7|) . Press and Schechter suppose 
that the eliminated volume consists of objects with mass between M and M + dM, 
corresponding to the idealisation that filtering the density contrast on any mass scale 
M cuts out precisely those objects with mass less than M while leaving unaffected 
objects with mass bigger than M. Ignoring the overdensity ~ (l + 5c) of the collapsed 
regions this implies that the number density dn of such objects is given by 
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Press and Schechter multiplied this formula by a factor 2, so that when integrated 
over all masses it would give the total mass density of the universe, rather than 
just the half corresponding to the regions of space where the linearly evolved density 
contrast is positive. They thus arrived at the estimate 



n(> M) ^ nps = / 

JM 



u e 



dM' 



(6.26) 



In this equation, R^ = Rf{M'), and similarly for < k"^ >' and z/'. The factor 2 inserted 
by Press and Schechter is not justified by their argument, because the linearly evolved 
density contrast has nothing to do with reality in the non-linear regime a{M) > 1. 
On the other hand, the neglected overdensity gives a factor ~ (1 + 5c) = 2 to 3. Thus 
the factor 2 goes in the right direction, and the Press-Schechter formula is reasonably 
well founded theoretically. 



From Eq. (p 



a{M) 
a{M') 



(6.27) 



The right hand side is independent of b, 6c and z, so it follows that rips, like ripeakj 
depends on these quantities through z/, which involves the combination Mc(l + z). 
In figure 15 is plotted the ratio of the two alternative estimates n{> M) = rips and 
n{> M) = npeak, for M = 10^° Mq and for M = IO^^Mq. One sees that the estimates 
agree to better than a factor 2 for z/ < 2. Presumably, this indicates that in this regime 
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the assumptions underlying the two estimates are compatible, in that increasing M 
by a small amount cuts out portions of the collapsed regions which have mass of order 
M and are centred on peaks with height of order z/(M). 

For large i/, the Press-Schechter estimate falls below npeak- This can be understood 
analytically, from the expression 

dripeak dnup _ dn^^ du dn^p d(/c^) dn^p d'j 
dM ~ dM ~ du dM d{k^) dM d-f dM ^ ' ^ 

The first term dominates for large u, leading to the ratio 



drips/dM ^ / 3 
dnpeak/dM [R}{k^), 



3/2 



i^--' (6.29) 



Apart from the factor 2, this is just the filter volume divided by the average volume 
of a collapsed region (Eq. ( |6.20|) ). 



Both estimates are too big when u is big. The estimate ?7.peak(— ''icoii) is too big 
because it counts many 'collapsed regions' which have mass less than the filtering 
mass and do not therefore contribute to n{> M). To obtain an accurate estimate, 
one would have to excise such regions before evaluating ripeak- The Press-Schechter 
estimate is too big because it assumes that all of the reduction in volume of the 
collapsed regions in going from M to M + 6M corresponds to structures with mass 
in this range, whereas some of the eliminated volume will belong to collapsed regions 
with mass less than M. To obtain an accurate estimate, one would again have to 
excise these regions before calculating the eliminated volume and hence obtaining 
the Press-Schechter estimate. Neither of these effects can be quantified analytically, 
because relevant statistical results are not known. The fact that the ratio Eq. ( |6.29|) 
is 2 times the ratio of the filtering volume to the average volume of a collapsed region 
presumably indicates that the regions which ought to be excised dominate both npeak 
and the reduction in volume leading to the Press-Schechter estimate. The implication 
is that even the latter is substantially too high at large u. 

In what follows we focus the estimate n{> M) = npeak(— n^p) because of its 
simplicity, but give some results also for Ups. 



6.6 Bias factors 

Before comparing the number densities with observation, we need one more result 
from the analysis of BBKS, namely an estimate of the galaxy and galaxy cluster 
bias factors at the present epoch. According to the CDM cosmogony such a bias 
factor will be present for any class of objects, if they form at an epoch when u is 
significantly bigger than 1 and therefore originate as exceptionally high peaks of the 
density contrast. In general the bias occurs partly through the excess clustering of 
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these peaks, and partly through additional non-linear clustering after the objects have 
formed. An estimate including both effects is 



a{M) 

where 



obj = 1 + T777T (6-30) 



u = u ^ — (6.31) 



1-72 

In these expressions, everything on the right hand side is to be evaluated at the epoch 
when the objects form. 



6.7 Comparison with observation 

The prediction n(> M) = ripcak is plotted against redshift in figure 16, for masses 
M = IO^^Mq, lO^^M©, 1O^°M0 and IO^Mq. For each case, three curves are given 
corresponding to the parameter choices {n,b} = {1, 1}, {1, 1.6} and {.7, 1.6}. Each 
curve ends at the epoch z^\{M), when the linear approach ceases to be valid. 

Let us consider first the case M = IO^^Mq, which corresponds to very large galaxy 
clusters. Since this mass corresponds to the normalisation scale Rf ~ lOh^^ Mpc 
there is little dependence on n in the range .7 < n < 1. There is however strong 
dependence on b, and several authors (Kaiser 1984; Bardeen et al 1986; Bardeen et al 
1987; Dalton et al 1992; Nichol et al 1992; Efstathiou, Bond and White 1992; Adams 
et al 1992) have claimed that 6 > 1.6 is preferred over 6 = 1. Let us first reproduce 
this result, then comment on its uncertainty. 

Since ^(IO^^Mq) = .8/6, we assume that the filtered density contrast on this scale 
is still evolving linearly at the present epoch. The quantities of interest are the number 
density ra(> M) and the ratio 6c/6 of the cluster bias factor to the galaxy bias factor. 
The prediction for n(> M) is 5.0 x IQ-^Mpc"^ if 6 = 1, and 1.0 x IQ-^ Mpc"^ if 
6 = 1.6. Following for instance BBKS and Bardeen et al (1987), we identify clusters 
in this mass range with Abell clusters of richness class > 1, and hence observed 
number density 7.5 x lO^^Mpc"^ (Bahcall & Soniera 1983). Thus, the prediction for 
n(> M) is about right if 6 = 1.6, but too big if 6 = 1. Coming to 6c/6, one finds 
from Eq. ( |6.30|) that if 6 = 1.6, then 6c = 4.0 giving 6c/6 = 2.5. If, on the other hand. 



6 = 1.0, then 6c = 1.5 which gives 6c/6 = 1.5. Recent estimates (Dalton et al 1992; 
Nichol et al 1992) give ^cc(r) = {ro/rf-'^^-^ with ro = (13 ± 5)/i~^Mpc. Dividing by 
the galaxy correlation function E,gg = {5h~^ Mpc/r)^'® gives therefore bc/bg = 2.4 ± .9 
as an observational estimate. Again, 6 = 1.6 is preferred. 

So far so good, but what about the uncertainty? Unfortunately it is big. First, 
suppose that we take 6c to be equal to 1.69, instead of the 1.33 used in the above 
estimates. For 6 = 1 (1.6) this multiplies n(> M) by a factor .56 (.16), and multiplies 
6c by a factor 1.4 (1.4). Second, suppose instead that we multiply M by a factor 
2, on the ground that there may be this amount of uncertainty in the observational 
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value of the masses of galaxy clusters (and remembering also that filtering the density 
contrast on mass scale M does not completely eliminate all structure with mass less 
than M). This multiplies n(> M) by a factor .30 (.09) and multiplies he by a factor 
1.6 (1.7). Third, do neither of these things but replace rzpeak by n^^ which is lower 
(though not low enough remember) for large v. This multiplies n(> M) by a factor .78 
(.35) without changing he- Each of these not unreasonable changes has a significant 
effect on the predictions. Finally, implement them all simultaneously. This multiplies 
n(> M) by a factor .03 (.0004), and multiplies he by a factor 2.3 (2.4). The effect is 
to make 6 = 1 strongly preferred over h = 1.6! 

Our conclusion is that one cannot yet reliably constrain h by considering the 
present number density and bias factor of rich clusters. On the other hand, these 
quantities are certainly very sensitive to h, and would constrain it well if the uncer- 
tainties could be removed. Similar remarks probably apply to earlier epochs, but as 
yet the data are too sparse to say much. 

For galaxies, the linear epoch ends before the present. The number density n{> 
M) at the end of the non-linear epoch cannot be estimated reliably, and it will also 
evolve with time because of merging and other non-linear phenomena. As a result, 
no reliable prediction is possible for the quantity n{> M) at the present epoch. 

However, in the biassed theory of galaxy formation, the formation of luminous 
galaxies is supposed to stop during the linear regime, at some epoch given roughly by 
Eq. ( |6.30| ). By demanding that this value reproduces the value of h used to normalise 
the amplitude one obtains the epoch of luminous galaxy formation, and hence the 
observed number density, as a function of h. This approach has been implemented 
for n = 1 by BBKS and by Bardeen et al (1987). Here we extend the calculation to 
n = .7, and comment on the uncertainty. 

The stars in figure 16b for h = 1.6 indicate the epoch of formation of luminous 
galaxies which is needed to reproduce this bias factor, according to Eq. ( [6.301 ). This 
epoch is 2; = 4.7 (2.8) for n = 1 (.7), and the corresponding value of n{> M) is 4.8 x 
lO-^Mpc"^ (7.9x lO-^Mpc"^). I gnoring merging etc., this comoving number density 
should be equal to the presently observed number density ng{> M) of luminous 
galaxies for M = lO^^M©. An observational estimate of %(> M) is the Schechter 
parametrisation, which taking the ratio of luminosity L to mass M to be independent 
of M is (Elhs et al 1988) 

In this formula 

0, = 1.56 X lO^^/i^Mpc (6.33) 

Integrating Mdug/dM over all masses gives = 1.6f2gai x IO^^Ii^^Mq, where Qg^i 
is the contribution to Q of luminous galaxies (including the dark halos) and we take 
^gai = For M = IO^^Mq, this gives an observational estimate ng{> M) = 
1.75 X 10"'^ Mpc""^, which is marked by an arrow on the y-axis of figure 16a. 
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According to this calculation, the predicted number density of luminous galaxies 
is somewhat too small for b — 1.6, and in fact one needs a value b ~ 1.2 to 1.3 to 
reproduce it. But now set M = 5 x lO^^M© in the theoretical calculation, on the 
ground that the observational value of M could be a factor 2 too high. With b = 1.6 
this gives a somewhat earlier epoch of formation z = 5.9 (3.6) for n = 1 (.7). The 
corresponding number densities are ng{> M) = 5.8 x 10^^ (1.2 x 10^'^), which look 
much more healthy, and the actual value of b needed for consistency is now around 1.4. 
Even without looking at any other sources of uncertainty (such as galaxy merging, 
which BBKS emphasise), it seems clear that it is very difficult to pin down b within 
the range 1 < 6 < 1.6 from these considerations. 

One conclusion, though, does seem fairly robust, which is that going from n = 1 
to n — .7 makes the redshift of bright galaxy formation later by around 2 units. As 
evidence accumulates from A'"-body simulations and from observation, this difference 
may eventually be able to pin down the required value of n. 

Finally, consider the results for M = IO^^Mq and IO^Mq, as shown in figures 16c 
and d. The observed number densities %(> M) according to the Schechter parametri- 
sation are again shown by arrows (though they become increasingly uncertain as the 
mass is reduced). Within the biassed galaxy formation theory one expects luminous 
galaxies with these masses to form before the epoch 2;ni(-^), but the observational 
bias factor for them is not known, and there is no reason why it should be equal to b 
which refers to bright galaxies. However, the epoch of formation would need to be a 
lot earlier than z-ai{M) to give agreement with the observed %(> M). It is presum- 
ably more reasonable to suppose that merging has taken place, reducing ng{> M) 
from its original value. Again, there does not seem to be a useful constraint on the 
parameters n and b. 

The quasar density 

So far we have focussed on observations at small redshift. One can also ask about 
observation at z >, 1. In particular, quasars have now been seen out to a redshift 
of about 5, and are the oldest observed objects in the universe. In order to produce 
the observed luminosity, some galaxies must have evolved to contain a sufficient con- 
centration of mass-energy. One then needs to estimate the minimum mass required 
and the number density, in order to discover if the CDM cosmogony can explain the 
observations. Such a comparison was made for n = 1 by Efstathiou and Rees (1988). 
Using both the ripeak and Press-Schechter prescriptions they calculated n{> M) for 
various masses as a function of redshift. Our results are far higher than theirs, but 
the difference can be explained if they chose b — 2.5 and defined b with respect 
to J3 normalisation instead of cr^ normalisation (they are not explicit about these 
choices, but our calculation essentially agrees with theirs if these changes are made 
in it). As in the case of large clusters at the present epoch, the extreme sensitivity to 
normalisation is caused by the fact that ly is substantially bigger than 1. 
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The number density of quasars assumed by Efstathiou and Rees was 10~^h^ 
Mpc~^, out to a redshift of around 4. Since then, observations have become more 
stringent, and to be consistent with present data one requires this number density 
out to a redshift of 5 (Martin Rees, private communication). This value also assumes 
that the quasar lifetimes are not too short, in which case one needs multiple genera- 
tions. One also needs to know the mass which is required to be evolving nonlinearly 
to harbour the quasar, and they estimated at the time that IO^^Mq was the smallest 
that one could safely get away with. However, recent work (Martin Rees, private 
communication) has suggested that only 10^°Mq is required, which corresponds to a 
large loosening of the constraint. As one sees from figures 16b or 16c, the upshot of 
all this is that enough quasars may form at high redshift, even if n = .7 and b = 1.6. 

Another cosmological requirement on high redshifts is the Gunn-Peterson con- 
straint on the amount of neutral hydrogen in the inter-galactic medium. In the CDM 
cosmogony this implies that some structure has formed before 2; = 5, in order to 
re-ionise the hydrogen (Schneider, Schmidt & Gunn 1989), but it is not clear how 
much or of what kind. Even with n = .7, structure with M < lO^^M© forms at the 
epoch 1 + z ^ 9/6 so a value for b as high as 1.6 need not be problematical. 



6.8 The virial velocity in a galaxy or cluster 

The linear approach has been pushed further by some authors, to try to predict the 
virial velocity v{M) of structures of mass M (Blumenthal et al 1984; Bardeen et al 
1987; Evrard 1989; Henry & Arnaud 1991; Evrard & Henry 1991; Adams et al 1992). 
The virial velocity of a gravitationally bound system is defined as the rms velocity 
of its constituents in the centre of mass frame. According to the virial theorem it is 
given by 

V = -— (6.34) 

Kg 

where M the mass of the system and Rg is its gravitational radius, defined by the 
requirement that its potential energy is —GM^/Rg. The idea is to relate Rg to the 
comoving size -Rcom of the object in the early universe, defined by 

M = (47r/3)i?,^„^po (6.35) 

where po = S-ffg/ (87rG') is the present mass density. To obtain such a relation one can 
use the spherical collapse model of Section 6.2, in the following manner. First, assume 
that the system has virialised, with no energy loss, before the epoch -Zform when in 
the absence of pressure forces the system would have contracted to a point. Next, 
assume that the system collapsed from an initial configuration at rest, with a density 
profile of the same shape as the profile after virialisation and without energy loss 
(dissipation). It then follows from the virial theorem, Eq. (|6.34 ), that Rg is equal to 



one half of the initial gravitational radius. Finally, set the initial gravitational radius 
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equal to the initial radius of the edge of the object, defined as the sphere containing 
mass M. Using the results already quoted in section ^2 this gives 




(6.36) 



(6.37) 
(6.38) 



Then Eqs. ( |6.34D and ( |6.35| ) give 
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This is the expression quoted by several of the authors mentioned above (Evrard 
1989; Henry & Arnaud 1991; Evrard & Henry 1991), and the others presumably used 
a similar expression though they are less explicit. Obviously, the assumptions leading 
to this relation are at best extremely crude approximations. It seems clear, in fact, 
that any unique relation between the mass, virial velocity and formation epoch of 
gravitationally bound systems can only be a rough approximation. 

In figure 16, the arrows on the horizontal axes show z^ormiv)-, calculated from 
Eq. ( |6.39| ), corresponding to estimates of the upper and lower observational limits on 
observational values of v for objects of mass M. These estimates of v are extremely 
crude. For galaxies, they correspond roughly to those given by Blumenthal et al (1984) 
without any attempt to update that analysis in the light of more recent observations. 
For clusters, the upper limit corresponds roughly to X-ray observations (Henry & 
Arnaud 1991); the observational lower limit of about 1000 km s~^ would in that case 
correspond to negative z, a value of about 1300 km s~^ corresponding to z = 0. 

The most direct way of utilising the relation Ziorm{M,v) is to set -Zform equal to 
Zn\{M). This gives the virial velocity v{M) of systems formed when the scale M goes 
nonlinear. Blumenthal et al (1984) identified such systems with luminous galaxies 
and with galaxy clusters, but according to the biassed galaxy formation theory that 
identification is wrong in the former case. Alternatively, the relation can be combined 
with a theoretical estimate of the number density dn{M, z) of objects with mass 
between M and M + dM, which form at epochs between z and z + dz, to calculate 
the number density n{z, > v) of objects with virial velocity bigger than v, which exist 
at redshift z. This is essentially the approach taken by the other authors mentioned 
above. 

We saw in Section 6.7 that even an estimate of n{z, > M) is subject to large un- 
certainties, with the result that it cannot be used to reliably constrain the parameters 
n and h within the range allowed by the QDOT and COBE data. An estimate of 
n(z, > v) inherits uncertainties of the same type, plus additional ones coming from 
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the use of the relation Zform{M,v). It is unclear how to quantifying the total uncer- 
tainty, but it seems likely that it will again prevent one from constraining n and b 
within the allowed range. 

6.9 Summary 

In this long section we have discussed galaxies and clusters, making the maximum 
use of linear theory despite the fact that these objects are gravitationally bound. 
Our discussion has extended earlier n = 1 calculations to the case .7 < n < 1. 
It has treated some points of principle not considered in earlier discussions, and in 
some cases has made for the first time a crude estimate of the uncertainty in the 
calculations. Despite, or perhaps because of, these improvements the conclusion is 
that they do not provide significant additional constraints on the pair of parameters 
n and b, beyond those already provided by the QDOT and COBE data. 

7 Discussion and Conclusions 

As small scale structure does not seem to provide any additional constraints over the 
large scale results of sections | and H, our final constraints are those illustrated in 
figure 11. From this, we conclude that for natural inflation one must have n > 0.70 
(which is a separate 2-sigma exclusion on two pieces of data). If one is willing to push 
all the observations to their limits and fold in extra uncertainties such as allowing 
h < 0.5, it is just conceivable that such a model allows a fit to the clustering data such 
as APM. The most important requirement for this is that the microwave fluctuation 
level should be at the top of its permitted range. One sees from the spread of the 
dashed lines on figure 11 that if the COBE result is exactly correct then the limit on 
n will immediately be pushed well above 0.8. 

For power-law and extended inflation, the limit is already n > 0.84 at high con- 
fidence. This clearly excludes the possibility that these models can fit the clustering 
data, and indeed we have already remarked (Liddle & Lyth 1992) that in fact this 
limit rules out simple forms of extended inflation, as they give a big bubble constraint 
(Liddle Sz Wands 1991) n < 0.75. (Changing from a cold dark matter cosmogony to 
a different choice also does not seem to salvage the situation (Liddle & Lyth 1992).) 
The strength of this constraint is due in the main to the large contribution of gravita- 
tional waves to the microwave anisotropies at low n; n = 0.84 is in fact coincidentally 
close to the break-even point where gravitational and scalar modes contribute the 
same amount. Note once again that the limit will strengthen dramatically if the mi- 
crowave result comes down. If the COBE result is exact, then that will push n well 
above 0.9. 

Assuming one takes the QDOT result seriously, then the bias is also strongly 
constrained with a 2-sigma maximum of around 1.6. For both types of inflationary 
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model, the lowest allowed n has this maximum bias. Hence the most optimistic power- 
law model, which may just fit the clustering data, appears to be a natural inflation 
style model with values of n ~ 0.70 and b ~ 1.6. In addition to an improvement on 
clustering issues, such a model also seems to deal adequately with the formation of 
structure (contrary to our more pessimistic assessment in LLS). 

A preference has been expressed for models with n in the range 0.8 to 0.9 (see 
eg Salopek 1992; Davis et al 1992; Lucchin, Matarrese & MoUerach 1992), because 
these models are claimed as allowing a higher bias than the n — 1 model. However, 
with the current error range on the COBE data our figure 11 shows little benefit in 
this at present, because even at n = 1 the COBE error bars allow quite substantial 
biasses at n = 1 (and remember this is without including that the spectrum from 
say a chaotic inflation model has a slight but not insigniflcant slope and also mildly 
signiflcant gravitational wave contributions). Even at n = 1 it appears that QDOT 
rather than COBE provides the upper limit to the bias. However, as soon as the upper 
limit to the fluctuations comes down, either by tightened error bars or independent 
measurement, this effect may become more relevant. 

To conclude, we have found strong constraints on the slope of the primeval spec- 
trum when generated by various inflationary models. Most extended inflation models 
appear to be ruled out completely. Power-law inflation is viable, but only for values 
of n > 0.84, too high to allow an explanation of the clustering data. Natural inflation 
(and related) models are the most promising candidates for generating useful power- 
law spectra — provided the true level of fluctuations is close to the top of the COBE 
range they seem marginally able to explain the excess large scale clustering as now 
seen in many optical surveys. In most aspects, such a model does at least as well 
as a bias one standard CDM model, with the advantage of a more plausible cluster 
abundance as well as helping with galaxy clustering statistics. We must note however 
that even without including corrections to the flat spectrum and a gravitational wave 
contribution, even standard CDM appears compatible with biasses up to 1.4 or 1.5. 
Such a model may however have difficulties with the clustering data. 

Final Note: As we were completing this paper, we received preprints by Cen et al 
(1992) which discusses various astrophysical consequences of power-law spectrum in- 
cluding A^-body results (but without considering gravitational waves), and by Adams 
et al (1992) discussing the astrophysics with particular emphasis on the natural infla- 
tion scenario. Where we cover the same ground as these papers, our results appear 
to be in good agreement. 
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Figure Captions 



Figure 1 

A selection of transfer functions, from Efstathiou (1990) (used in this paper), Bardeen 
et al (1986), Davis et al (1985) and Holtzmann (1989), given as appropriate to CDM 
universes with Q = 1, /i = 0.5, low baryonic content and three neutrino species. They 
feature 4, 5, 3 and 4 fitting parameters respectively. 

Figure 2 

The normahsation bSnik — 1 Mpc~^) as a function of n. The solid line is from 
the CTg normalisation (used in this paper), while the dashed line shows the rival J3 
normalisation. 

Figure 3 

The present day spectra, as calculated in the linear approximation, for a selection of 
values of n. One sees the additional large scale power and the deficit on short scales 
when one compares n < 1 to the standard CDM spectrum. Note k is in Mpc~^, 
without a factor of h. 

Figure 4 

The dispersion ba{M) as a function of mass, for n = 1 and n — 0.7, and with both 
choices of filter. The top hat filtered spectra are unity at M ~ lO^^M©, as required 
by the normalisation. The gaussian filtering gives significantly smaller answers than 
does the top hat, as its smearing gives a higher contribution to the larger scales. 

Figure 5 

The triangles indicate the prediction for the mean quadrupole as a function of n for 
power-law inflation, and the stars for natural inflation. The dotted lines plot the 
COBE observation for two choices of bias (note the bias implicit in the |/-axis). The 
CDM prediction for the mean at bias 1 is very close to the COBE result. The vertical 
bars on the starred data indicate the spread of the pdfior the quadrupole (for clarity 
the bars have been omitted for the triangles — they are the same size); 95% of the 
pdf is above the bottom of the bars, while 95% is below the top. The xl distribution 
is not symmetric, so the bars are skewed to higher values (somewhat concealed by 
the log plot). A value of n is allowed at 95% exclusion if the observations cut through 
its vertical bar. ModeUing the observational errors (see text) gives an even looser 
criterion. 

Figure 6 

The predictions for b ^ as defined in text, along with the COBE limits (at 1- 

sigma) for bias 1 and 2. The upper line represents the power-law inflation predictions, 
the lower those for natural inflation. As for the quadrupole, the n = 1 prediction at 
bias 1 is very close to the COBE result. 
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Figure 7 

The predicted angular correlation functions for a choice of n are plotted alongside the 
observational data from the APM survey (Maddox et al 1990). With the anticipated 
residual systematics, values of n between about 0.3 and 0.6 provide reasonable fits, 
while the standard CDM curve falls well below the data. 

Figure 8 

A comparison of models with different n against a standard CDM model but with 
Vth = 0.2. 

Figure 9 

The predicted rms velocity flows, when smoothed with a top hat of radius i?/, for 
different choices of n. 

Figure 10 

The predicted rms velocity flows in a configuration mimicking the POTENT observa- 
tional data. The velocity field is first smoothed with a 12h~^ Mpc gaussian, reducing 
the short scale power, and then smoothed with top hat filters of radius i?/, giving 
predictions smaller than in figure 9. The solid lines indicate the predictions for n = 1 
and n = 0.7. The stars indicate the POTENT observational data at bias 1 (read 
from figure 4 of Dekel (1991)), and the triangles the same at bias 1.6. The error bars 
on the data (the last ones just overlap) are 1-sigma. Finally, we emphasise that the 
theoretical curves are averages over all observer points, whereas the observations are 
a single realisation, with correlated errors due to long wavelength domination of bulk 
flows. 

Figure 11 

The limits on the bias 6 as a function of n. The dotted line shows an upper limit on n 
to fit the APM data. All other lines are 2-sigma limits. The dot-dashed lines indicate 
the range allowed by QDOT, and apply to all models. The solid and dashed lines 
indicate the limits from COBE for power-law (and extended) inflation and natural 
inflation respectively. For power-law inflation, one has n > 0.84 and b < 1.5, which 
rules out extended inflation which needs n < 0.75 to satisfy the bubbles constraint. 
For natural inflation, one has n > 0.70 and b < 1.6. This is well outside our APM fit 
range (though it is close to the range allowed by Efstathiou, Bond & White (1992)). 
If one does not attempt to fit optical galaxy clustering, then the allowed region is 
that marked with the star, bounded by the microwave limits appropriate to a given 
scenario. 

Figure 12 

The asphericity parameter x"^ plotted against v. 
Figure 13 

The ratio 2n-a^/n-^ plotted against u. 
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Figure I4 

The ratio l^peak/^/ plotted against u. 

Figure 15 

The ratio nup/np_s plotted versus v. The ratio is not very sensitive to n, and only 
the case n = 1 is plotted. 

Figure 16 

The number density n{> M) as a function of redshift. It is calculated theoretically by 
equating it with the number density of peaks of the linearly evolved density contrast, 
filtered on the mass scale M. Figures 16a-d refer respectively to M = lO^^Mo, 
IQ-'^^AfQ 10^*^ M0 and lO^M©. For each case there are three curves. They correspond 
to the three choices {n, h} = {.7, 1.6}, {1, 1.6} and {1, 1} which mark the corners of the 
more or less triangular region allowed by the QDOT and COBE data (figure 11). The 
arrows on the vertical axes give the observed galaxy number densities at the present 
epoch. The arrows on the horizontal axes indicate very roughly the 'observed' range 
of formation epochs, deduced from the indicated range of virial velocities. Each curve 
ends at the epoch when a{M) = 1, signalling the end of linear evolution. According 
to the theory of biassed galaxy formation, luminous galaxies form significantly before 
that epoch. In the case of bright galaxies, the epoch can be calculated by demanding 
that the resulting bias factor is equal to b, and it is indicated by a star in figure 16 
b. The theoretical and observational uncertainties in are discussed in the text. 



51 



